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Education and research in 


ENSINO E INVESTIGAÇÃO EM MECÂNICA NOS ESTA- 
DOS UNIDOS 


O artigo dá informação sobre a organização do 
ensino e investigação no dominio da Mecânica nos Es- 
tados Unidos. Foi apresentado pelo autor durante a 
última sessão, dedicada ao Ensino e à Investigação em 
Mecânica, do 1.º Congresso Português de Mecânica 
Teórica e Aplicada [Dezembro de 1974). 


[Um supposed to talk about the organization of Edu- 
cation and Research in Mechanics in the USA. However, 
there is no such organization. Unlike most European 
countries, there is no Central, Federal, Ministry of Edu- 
cation that has control of eithe: education or research. 
In fact, there are only three federal schools in the USA: 
the schools that train the army officers, the navy officers 
and the aviators, and there is no Mechanics ;esea:ch or 
higher Mechanics education conducted in those schools. 

There are about 3.500 institutions of higher edu- 
cation in USA. Most of them are small, privately 
endowed schools. In addition, each of the 5O states 
which make up the United States has at least one and 
sometimes many more institutions of higher education, 
For example, New York State, in which New Yo:k 
City is situated, has about 5O institutions of higher 
education of one kind o: another, and California 
has probably as many. Some of the states of smaller 
population have fewer but there is least one ins- 
titution of higher education and each state maintains 
a public institution that is suppo:ted mainly by the taxes 
from the State. 

In further addition, there are many private Univer- 
sities and private Institutes. Most of the names are pro- 
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Professor 


EDUCATION AND RESEARCH IN MECHANICS IN THE 
U.S.A. 


Information is given about the organization of Edu- 
cation and Research in the field of Mechanics in the 
U.S.A. The paper was presented by the author during 
the last session, devoted to Education and Research in 
Mechanics, of the First Portuguese Congress on 
Theoretical and Applied Mechanics [December 1974). 


bably familiar to you: Ha:vard, Yale, Princeton, Colum- 
bia, Stanford; these a:e all privatelly endowed schools 
which ceceive essentially no support either from the 
states or from the federal government. Then there are 
a few Institutes of Technology, whose names are pro- 
bably also familiar to most of you: «California Institute 
of Technology», «Massachussetts Institute of Techno- 
logyh. There are also two othe: prominent ones: «Brook- 
lyn Polytechnical Institute», «Stevens Polytechnical Ins- 
titute», 


These Institutes and the Engineering Schools in the 
Universities, either private or public Universities, do 
receivo funds indirectly fiom the Federal Government, 


To those institutions that do research (for example, 
in Mechanics, there ae about 100 institutions in which 
research in Mechanics is conducted) federal funds are 
supplied through individual Professors. The individual 
professor will apply to one of perhaps 5O different go- 
venment agencies for a contract to support his re- 
search; and, if he is successful, his university gets 
certain funds. Mainly those funds go to supply fellowhips 
for students to study for advanced degrees. Part of those 
funds, sometimes, goes to pay a part of the professor's 


(*) Paper presented at the 1 st. Portuguese National Congress on Theoretical and Applied Mechanics, Lisbon, December 1974, 


(**) Professor, Columbia University (New York). 
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salary and also a part of the funds might go towards 
buying the equipment necessary for laboratory work. 

Now, of the 100 institutes or universities or schools 
that have departments in which mechanics is taught on 
an advanced scale, there is complete independence. In 
each Department of Mechanics or Department of Civil 
Enginee.ing in which mechanics may be taught, or De- 
partment of Mechanical Engineering in which mechanics 
may be taught, each Head of Mechanics is King in his 
own Empire; but he is not in charge, because there are 
p:incipalities in each one of these 109 Kingdoms. There 
is a Prince of Solid Mechanics, a Prince of Aerodyna- 
nics and more. Each one of whom usually has complete 
control of what he teaches and what research he does; 
and, so, you will find a 100, times maybe seven, or 
700 principalities, each independent. The courses of 
instruction will vary quite widely from one extreme to 
another. The type of research is the choice of the Prince; 
and in a Civil Engineering Department you might find 
one Prince who is working in electromagnetic theory 
because that is what he happens to be interested in. 
There is ve:y great freedom and so there is a very wide 
variety of both instruction and research over the whole 
country. 

Now, there is no Central Administration. However, 
a few yea:s ago, as result of a combination of circums- 
tances, (President Nixon did not like scientists and so 
abolished «The President's Science Advisory Board») 
there came a tremendous reduction in the government 
funds available to individual professors to conduct re- 
search. (These funds also help the school, because they 
pay part of the professor's salaries. There is also what 
we call «overheady which goes to the English Depart- 
ment or the History Department or to help pay other 
depa;tment or university expenses). As result of this re- 
duction in funds, the Engineering Schools looked 
around: what shall we drop? Well, the easiest thing to 
drop is Mechanics, because nobody quite knows what 
is Mechanics: is it Civil Engineering, is it Mechanical 
Engineecing, is it something else?. And, as result of 
that, the 100, actually the 101, Chairman of Departments 
of Mechanics got together and formed a protective orga- 
nization. lt is called «The Association of Chairman of 
Deparitments of Mechanics» and this is really the first 
gathe.ing together of Mechanicists, from the educational 
view, in the USA. There are 101 members and | have 
tu-ned over to your chairman the list of names, 
addresses end telephone numbers of each one of them; 
so, it Is easy to get in touch with them, 


This Association publishes a quartely News Letter 
wich gives news of what is going on Mechanics Depart- 
ments throughout the country and they also publish 
Proceedings of meetings on special topics. For example, 
there was a meeting on the topic of the place of Me- 
chanics in the engineering curriculum and Il give a 
copy of that to your Chairman. That Association was 
organized by Bruno Boley who is a great mover, a very 
activa member of our mechanices community. He also 
started, at the seme time, an organization called «The 
American Academy of Mechanics» and this has about 
500 members and it publishes a monthly jou-nal, some 


copies of which | have also turned over to you chairman. 
That monthly journal is a littla more technical than the 
Association of Chairman journal is. 


Now, there is one organization which presumably 
ovesees the courses of instruction in Engineering 
Schools and that is called «The Engineer's Council for 
Professional Development», that is a council composed 
of the five representatives from the five major Engi- 
neering Societies of the country: «The American Society 
of Civil Engineers», «The Amesican Society of Mechani- 
cal Engineers», «The American Society of Chemical En- 
gineers», «The American Institute of Elect.ical and 
Electronics Engineers» and «The American Institute of 
Mining and Metalurgical Engineers». They have appoin- 
ted committees which are supposed to accredit Enginee- 
ring Schools. They send around a small commiitee to 
Bach engineering school and this committes submits a 
report saying this engineering school has, for example, 
an acceptable Department of Civil Engineering; or they 
might say, well, in this Department of Civil Engineering 
they a:ie very strong in Mechanics of Solids but they 
ara very weak in Mechanics of Fluids. This is usually 
what occurs, for example, if a representative on the 
inspection committee happens to be a man in Fluid Me- 
chanics. The standards of ECPD (Engineers Council fo; 
Professional Development) are not ve:y high; it is very 
rare that the Department would be not acc:edited; but 
that is a commitee formed from five private societies 
and there is no connection with the government 
whatsoever. 

There is also a Society fo: Engineering Education 
and this has a Mechanics Division but it is concerned 
mainly with problems of administration: little with 
education. 

So much for the «organization» of education in 
mechanics in the USA. Thera is no organization, it is 
completely haphazard and, personally, | think that is a 
good thing because it makes for very broad coverage. 
There is no limitation on what might occur; and this 
gives ample room for expansion in all directions. 


And, now, as for o-ganizations concerned with 
research, these are essentially the technical and scien- 
tific Societies. As far as Mechanics in concerned, the 
largest and the most active is the Applied Mechanics 
Division of the American Society of Mechanical Enai- 
neers. There are about 11.000 members of that Division. 
The reason that Division is so very strong and p.edo- 
minant, is that it was the first foothold that Timoshenko 
gained when he came to the USA. There was then little 
Applied Mechanics in the USA and, as he was a Civil 
Engineer, he tried to start a Division of Mechanics in 
tho American Society of Civil Engineers, but they 
would have none of it. And so he went to the Mecha- 
nical Engineers and, as result, there is now the most 
active division or society that has to do with research in 
Mechanics. The Engineering Mechanics Division of the 
American Socieiy of Civil Engineers came along some 
dozen years late: and that is now the next most active 
having to do with research in Mechanics. 

Both those Societies publish journals. The ASME 
(Mechanical Engineering Society) publishes the Journal 
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of Appliad Mechanics, Applied Mechanics Reviews, and 
Applied Mechanics News; and the Civil Engineering So- 
ciety publishes the Proceadings of the Engineering Me- 
chanics Division, and also a Mechanics News Bulletin. 

Next, in activity, is the Society for Experimental 
Stess Analysis; it is a relatively new Society: it's 
only about 30 years old. They a very active group 
in both Experimental and Theoretical Mechanics. They 
publish a monthly journal «Experimental Mechanics» 
and an annual «Proceedings». 

Then comes the «American Institute of Aeronautics 
and Astronautics», which is quite active in meetings 
and publications in those aspects of Mechanics which 
apply to Aeronautics and àAstronautics, in particular 
aerodynamics and airc:aft structures, and these papers 
appea: in the AIAA Journal with which many of vou 
are familiar. Also the American Institute of Aeronautics 
and àAstronautics co-sponsors one meeting each year 
with Mechanical Engineering Society on Structures and 
Materials problems. 

Then there is a small Society of Engineering 
Science which holds one annual meeting and whose 
founder edits the monthly «International Journal of En- 
gineecing Science. The same publisher, Pergamon Press 
publishes the International Journal of Solids and Structu- 
res, Computers and Structures, Engineering Fracture Me- 
chanics and Mechanics Research Communications. 


And, finally, there is the Acoustical Society of 
America, whose main program is acoustics. The papers 
presented at their meetings are published in the Journal 
of the Acoustical Society. À substantial portion has to 
do with Mechanics: vibration problems, wave propaga- 
tion problems, ultrasonics problems. This is one of the 
main journals in which mechanics research appears. 

Then, there are two more organizations, they are 
realy just committees: «The Midwest Mechanics Confe- 
rence» and «The Southeast Mechanics Conference», 
each one of these holds annual meeting on mechanics 
and publishes Proceedings. 

Finally, there is the US National Committes on 
Theoretical and Applied Mechanics, that is the Com- 
mittee analogous to the one which has organized this 
meeting. This US Committee was organized about 24 
years ago they sponsor a National Congress every four 
four years. The Seventh Congress was held just this past 
summer. That committee is composed of one represen- 
tative from each of 10 societies: the American Society 
of Mechanical Engineers, The Amesican Society of Civil 
Engineers, the Society for Experimental Stress Analysis, 
The American Institute of Aeronautics and Astronautics, 
The American Mathematical Society, The American 
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Institute of Chemical Engineers, The American Physical 
society for Testing Materials, The Society of Rheology 
and The Society for Industrial and Applied Mathematics. 
Each one for those 10 societies has a representative on 
the US National Committes for Theoretical and Applied 
Mechanics and that's the Committee that is represented 
in IUTAM, the International Union. 

The US National Committee, mainly represents the 
US in IUTAM and chooses, among applicants, the Uni- 
versity which will hold the National Congress; and once 
it does make that choice its duty is finished. For exam- 
ple the last Congress was held at the University of Co- 
lorado and the University of Colorado is the one that 
o:ganizes the meeting, which sends out the announce- 
ments and does eve:ything in connection with the ope- 
ration of the meeting. The US National Committee 
supplies aproval, At the Seventh Congress that took 
place this past summer there were about 180 papers 
presented. About 70% of those papers were from univer- 
sity authors, about 15%, from industrial laboratories and 
about 15% from Government laboratories. So that you 
sea that, roughly speaking, about 70% of research in 
Mechanics is done in the Universities, about 15% indus- 
trial laboratories (mostly in large industries like the Bell 
Telephone Laboratories or IBM), and about 15% in 
government laboratories. 

The host University, at one time, published the 
Procedings; but, with inflation, it got to be too expensivea 
and, at the last Congress, instead of having published 
Proceedings, the Congress Committee, The Host Com- 
mittee, got together with forty-four co-operating Jour- 
nals: here is the list of them. These Journals agreed to 
publish the papers which are acceptable to them; so that, 
to get a broad view of what goes on in Mechanics re- 
search in the USA it is no longer possible to look at 
the Proceedings. What you must look at is the program 
of the meeting and then write to the author, whose 
address is given, and find out which Journal his paper 
is published in, and then look at the Journal. It is a 
little Indirect but, you see, there is no government 
support to these Congresses, The entire expense must 
either be paid by the host university or sometimes 
there were collected some small funds from industry 
o:, perhaps, some small support from a government 
agency. 

Weill, that is an overall view of the disorganization 
of Mechanics in the USA. It is hardly necessary to say 
that, despite the absence of organization, or perhaps 
because ot it, there are many centers of research at 
which work of high quality and originality is in progress. 
Thank you. 
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Some applications of functional analysis in the 


mathematical theory of structures O 


RESUMO 


As analogias formais entre modelos estruturais, 
assim como o problema da geração de modelos a pertir 
de outros modelos, caem dentro do âmbito da Teoria 
Matemática das Estruturas a qual pode ser portanto for- 
mulada em três partes: um modelo genérico, regras para 
gerar modelos a partir de outros modelos, uma justifica- 
ção pera tais regras. O presente trabalho considera es- 
sas diferentes partes, recorrendo à Análise Funcional 
como o quadro matemático adequado à formulação da 
teoria. 


ft. A modern view of the Theory of Struclures 


Let us start by explaining what the Theory of 
Structures really means to the author as a part of Solid 
Mechanics (seo [1], [2], [3]). 


Solid Mechanics comprehends different models 
conceived for the equilibrium and deformation analysis 
of solids, namely three-dimensional models, two-dimen- 
sional models, one-dimensional models and discrete 
models. 

The methods of solution of particular problems 
within the frame of each model do not fall within the 
scopo of the Theory of Structures but of other theories, 
like the Theory of Elasticity, the Theory of Shells, the 
Theory of Roads, and so on. 


The formal analogies between models, togethei 
with the generation of models from other models, fall 
howeve: within the scope of what the author calls the 
Mathematical Theory of Structures, which therefore may 
bo formulated as consisting of three parts: 


— à generic model, 
— rules for generating models from other models, 


— a justification for such rules, 


E. R. Arantes Oliveira (1) 
Technical University of Lisbon 


SUMMARY 


The formal analogies between structural models, 
together with the generation ot models from other mo- 
deis, fall within the scope of the Mathematical Theory 
of Structures which therefore may be formulated as 
consisting of three parts: a generic model, rules for 
generating models from other models, a justification for 
those rules. The present paper considers all such paris 
and resorts to Functional Analysis as an adequate ma- 
thematical frame for formulating the theory. 


Only elastic structures under static equilibrium will 
be considered in the present paper. Physical linearity 
will not be required. 


2. The generic model 


The generic model consists in three groups of 
equations — force-st-ess, strain-displacement and stress- 
-strain equations — involving four kinds of magnitudes 
—— stresses, strains, forces and displacements. Such 
equations are supplemented by force and displacement 
boundary conditions (external and eventuaily internal) 
and supposed such that the work principle holds. 

A couple of a stress and a strain field related by the 
st.ess-stiain equations form a structural field. A structu- 
ral field is called a compatible field with respect to a 
system of incompatibilities (a set of prescribed initial 
strains and displacement boundary conditions), or is 
said to compatibilize such system of incompatibilities, 
if it satisfies the strain-displacement conditions and the 
displacement boundary conditions. À structural field is 
called an equilibrated field with respect to a system of 
external forces (a set of prescribed body forces and 
force boundary conditions), or is said to equilibrate 
such system of external forces, if it satisfies the force- 
-stress equations and the force boundary conditions. A 


(*) Paper presented at the «Joint IUTAM/IMU Symposium on Applications of Methods of Functional Analysis to Problems of 


Mechanics», Marseilles, September 1975. 


(1) Professor Catedrático de Engenharia Civil, Investigador do CEMUL 
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field which is simultaneously a compatible and an equi- 
librated one is an exact solution (to the structural equa- 
tions and boundary conditions) with respect to given 
systems of incompatibilities and external forces. 

Let X be the set of all the structural fields associa- 
ted to a given elastic structure. 

Set X can be made a vector space by defining 
the operations of addition and multiplication by a scalar. 
Such definitions can be made in different ways among 
which the most natural ones are schematized below: 


x, corresponds to e, 


Xu » + 6, 
mo “Foxs » » ej+te, 
x, corresponds to s, 

Xa » +» 5, 
do ) » S,+5, 


(2.1) 


whe:e e and s represent the strain and stress vectors. 


X can still be made a Banach space by associating 
a norm to each of its elements. Such norm can also ba 
defined in several ways as, for instance, 


E ES 
par 4/ [ ecnda p= [s.sas 


Dr 


á 1 
(2.2) 


where A denotes tha domain corresponding to the stru- 
cture. 


| compatibilizel 
Tha set of all the fields in X wich o 
equilibrate 


| incompatibilities E 
a given system of | | is called an 
external forces | 


É 
| isocompatible | 


E | subset of X. 
isoequilibrated | 


" isocompatible 


The set of all the | subsets of X 


isoequilibrated | 


is assumed also a Banach space and denoted by 


J 
EI 


corresponds to each 


) | | incompatibilities 
7 - À system of | j 


external forces 


corresponds thus to each element of 


E 
Às a unique element of | f | 


element of X, a function + assumed continuous, 


which associates to each element of X tha corresponding 


with domain X and range , Can be considered 


element of | - Wa writa therefore 


[É 


LbJ=i | et)=E (2.3) 


Ie) 
E cf 


where 


Equations (2.3) left and right are respectively 
called compatibility and equilibrium equations. 


Assuming that the intersection of each isocompa- 
tible and each isoequilibrated subset of X contains one 
and no more than one element, a one-to-one correspon- 
dence can be established between the elements of X 
and the elements of the cartesian product |» £. Such 


cartesian product will be called the space A of the 
external actions while X may be called the space of 
responses. 


Although | and « have no inverse, a function X can 
still be considered with domain A and range X,which 
associates to each pair (Il, E) of an isocompatible and 
an isoequilibrated subset the corresponding unique 
intersection, x. We write 


x = x (1,E) 
(2.4) 


and assume x a continuous function. 


With the help of the work principle, the total po- 
tential and complementary energy theorems can be pro- 


| | potential 
ved. The total energy theorem sta- 


complementary 


| potentl- 
tes that the exact solution makes the total ê 
comple- 
mentary | 
» | energy stationary on the set of the 
' compatible | 


| E fields. Such theorems become minimum 
“ equilibrated 


theorems if stability is admited, 


complementary | 
The total | | energy being a con- 


potential 


Te be) 


tinuous functional, ú 
É; (x) 


E 
, on X (index | ; 


external fo | 
refers to the system of | ERR | 


AR which the 
| imcompatibilities | 
functional corrsponds to), the minimum total | Po- 


| com- 
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tential 
energy theorem may be enounced 
plementary 
E ce oà : | T E (x) om | 
by stating that x (|, E) minimizes ; 
Tylx) on E 


The distance between two structural fields, |. e., 
two elements of X, could of course be defined as the 
no:m (defined as above) of their difference, As it will 
become clear later, however, a more convenient defi- 
nition consists in making 


dês =1/, | ós. de dá 


(2.5) 


for two elements such that the norm of their difference 
is small, and then defining the distance between two 
arbitrary points x, and x, as the smallest obtained by 
adding the infinitesimal distances along all possible 
continuous (*) pathes in X between x, and x,, À. e. by 
making 


I (2.6) 


where ds denotes the distance, defined through (2.5), 
between two points very near each other. 


Às it may be proved that 


sm. | a 


6º Tp= | 5s. seda on l 6º Ty= [5s. seda on E 


A 


A 


des 


(2.7) 


at x=x (|, E) and as, at the same point, Te = 


=3T,/=0on , there follows, considering (2.5) 
that, if º | and hl is very small, 
hSE, | E 
o? (x + hxj= dº (x + h, x) = 


where E, and ls denote the zero elements of É and J. 


A more general scheme may be introduced [4] 
which covers any situation in which an extremum prin- 
ciple is known to exist. Such scheme simply considers 
space X and a family “&2 of continuous functionals q 
which are assumed to admit a proper minimizer s on 
each subset C belonging to a certain class of subsets 


of X, called constrained subsets of X.s is assumed a 
critical point of q on C,i. e., at point s, grad qy vanishes 
on €C. 

The constrained subsets are assumed homeomor- 
phic to a certain linear subspace of X and their union is 
assumed to coincide with X. 

The set of all the minimizers corresponding to all 
the constrained subsets of X is called a minimizing 
subset of X. Each minimizing subset corresponds to a 
certain functional q «1, and the union of the minimi- 
zing subsets corresponding to all the functionals of 
family & is also assumed to coincide with X. 

The intersection of each constrained and each 
minimizing subset of X is assumed to contain one and 
no more than one element. 
| constrained 


minimizing | 
| 


Elements belonging to the same 


i strain é 
inqconatra indo in. 
isominimizing 


subset of X will be called 
Reference to such general scheme is made, 
for sake of commodity, along the next Sections. 


3. The generation rules 


In what concerns the rules for generating models 
from other models, their importance in the Theory of 
Structures can be realized if it is considered that dis- 
cretization, i. e., the generation of discrete models from 
continuous models, is the typical method of the structu- 
ral analyst. As, on the other hand, the generation of a 
discreto model is quite analogous to, for instance, the 
generation of a two-dimensional model from a three- 
-dimensional one, there is no reason why the study of 
discretization should not be included within a general 
theory of the generation of models, which comprehends 
the description of general generation rules and their jus- 
tification with the help of convergence analyses. 

The procedures for generating new models in the 
Theory of Structures are essentially two dual ones: the 
potential and the complementary energy methods [1, 
2, 31. 


| potential 
The energy method starts by 


complementary | 


ia strains and displacements 
defining the generated 


stresses and tractions 


| | strain- 
in terms of the generating ones. Generated | i 
orce- 


displacement 
-Stre55 


displacement 


equations and 
"force 


boundary conditions are simultaneously introduced which 
must satisfy the condition of being exact in the frame 
of the generating model. Such equations introduce a 


(*) Continuity gained of course a meaning as 5s00n as the norm was defined. 
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correspondence between the generated and the gene- 
rating fields which may be represented by the equation 


Y =B,x 
(3.1) 


where B, is a linear operator whose domain is the space 
X of the generating fields and whose range is the space 
Y of the generated ones. 


As the generating and generated models are analo- 
gous to each other, constrained and minimizing subsets 
can be considered in Y as well as in X. The fact that 


compatibility 


the generated ndo 
equilibrium 


equations are exact 


in the generating model implies that the B,-images of 
isoconstrained elements in X are isoconstrained in Y. 

The work principle, or some variational principle 
resulting from the work principle, is subsequently used 


equilibrium 


for deriving the generated | id 
compatibility 


equations. 


An assumption is then introduced, like the Kir- 
chhoff's assumption of the theory of shells, or the 
Bernouilli's assumption of the theory of rods, or “he 
definition of the fields allowed within the elements in 
the finite element technique, which establishes a corres- 
pondence between the elements of Y and elements of 
X. A new linear operator, B,. is thus introduced with 
domain Y and range X', X' being a subspace of X (the 
space of the allowed fields), and we write 


x = B, ly) 
(3.2] 
potential | 
The energy method postulates 
| complementary | 
potential 
the invariance of the total energy. 
| complementary 
This means that 
TIB.(y]=V(y) | THB, (y)] = Vº (y) (3.3) 
Tand V PR; 
where | denote the total potential 
| T*andV* | complementa:y | 


energies in the generating and in the generated models, 
and permits to express the generalized elastic coeffi- 
cients and external force defining magnitudes in terms 
of tha generating ones [11]. 
According to what was written above and namely 
to equation (3.3), the linear operator 
B=8B,8B, (3.4) 
with domain X and range X', which we call the inter- 
polation, operator, enjoys the following properties: 


ij the B-image of any element belonging to X' 
coincides with the element itself, i. e., 


Bi(x)=x fx E A 


(3.5) 


ii) Constrained and minimizing subsets meeting 
the same requirements as those in X can be 
defined in X', with respect to the same family 
of functionals, b, 


ii) The B-images of isoconstrained elemenis in 
X are isoconstrained in X', although not ne- 
cessarily in X. 

lf any two isoconstrained elements in X' are iso- 
constrained in X', opecator B is said to be conforming. 

The constrained subset C' which contains tha B- 
-images of the elements of a given constrained subset 
Cof X is said to correspond to C. À given minimizing 
subset D' of X' is said to correspond to a certain mini- 
mizing subset D of X if they both cor;espond to the 
same functional q. 

A second operator, À, also with domain X and 
range X', can still be considered which makes the 
intersection of each constrained and each minimizing 
subset of X correspond to the intersection of the 
corresponding constrained and minimizing subsets of X'. 

Operator A, called the approximation opecator, is 
assumed bounded and continuous. The A-image x, 
of an element x=X is called the approximation of x in 
K'. 

We remark that the A-images of any two isocons- 
trained elements in X are isoconstrained elements in 
X', and that the A-images of any two isominimizing 
elements in X are isominimizing elements in X'. 

Consideration of operators À and B permits to 
discuss the conversion of a given variational problem 
into a new one, eventually simple; than the first and 
capabla of p-oviding an approximate solution to it. 

H the interpolation operator B is chosen in such 
a way that X' is N-dimensional, then, the variational 
problem becomes particularly easy to solve, as functio- 
nal q becomes a function of N variables, and the 
element which minimizes the functional can be deter- 
mined by just solving the system of N equations to N 
unknowns which results from equating to zero the deri- 
vatives of the function with respect to each variable. 
The problem is said to have been discretized. Of 
course, the equations will be linear only if functional 
y is a quadratic function of the variables. 

If operator B is conforming and C'- C, this dis- 
cretization procedure is nothing else than the classical 
Ritz method. In the case of the finite element method, 
however, operator B is generally non-conforming. 


4. dJustification by convergence 


The evolution of the modern Theory of Structures 
was deeply influenced by the finite element method, 
end it was namely in connexion with the finite element 
method that the role of convergence in the theory was 
fully appreciated, 

The importance of convergence (to the exact solu- 
tion) was truly realized by the finite element analysts 
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not ea:lier than 1965 [5]. Before 1965, indeed, too 
much stress was put upon confo:rmity and monotonic 
convergence, so that the search for conforming elements 
was then a popular topic of research. 


Later, if was realized that conformity is not strictly 
necessary and that it can even be inconvenient. The 
main theoretical reason for a given type of element 
baing accepted is thus really its capacity for generating 
sequences of approximate solutions tending to the exact 
solutions of a sufficiently wide class of problems. 


lt cannot be forgotten however that, if confermity 
is violated, the finite element method, regarded as a 
variational technique or, more precisely, as a technique 
which provides an approximate solution to a variational 
problem consisting in minimizing a certain functional q 
on a certain space C, by just minimizing q on a space C' 
not generally contained in C, becomes distinct from 
Ritz's method, which supposes that C' is contained in 
C. The convergence theo:ems of the last [6] could not 
the-efore be used with the fo:mer and a new convergence 
theory had to be built. 

The generation of discrete models by the finite 
element technique being thus justifigd by convergence, 
it is natural that convergence considerations be made 
clso for justifying the generation methods even if thêy 
aro used fo: the generation of continuous models Like 
the theory of shells. And this explains how the Theory 
of Structures can give answers to problem, like the cone 
of the foundations of the theory of shells, which have 
been seen, until quite recently, exclusively from other 
points of view (see [7]). 

The statement that the generation methods are 
justifigd by convergence requires some explanation, 
however, because, although it is clear what convergence 
means in the case of the finite element method, it 
seems less clear what it means in the case of the gene- 
ration of the theory of shells, for instance. 

Let us consider then the case of a shell and assume 
that its thickness is made smaller and smaller. What 
the convergence analysis is required to prove is that 
the two-dimensional (generated) solution becomes more 
and more near the coresponding three-dimensional 
(exact) solution as the shell becomes thinner and 
thinner, provided the relative values of the bending and 
membrane elastic coefficients (generated elastic coeffi- 
cients) are not changed when the thickness tends to 
zero, 

If such coefficients really do not change, then, 
tha two-dimensional (generated) solution does not change 
also and what must be appreciated is the convergence 
of a sequence of generating (th: ee-dimensional) solutions 
towards a generated one, and not, like in the finite 
element case, the convergence of a sequence of gene- 
rated (app-roximate) solutions towards a generating 
(exact) one. 

There remains to remark that the condition of 
unchangeability of the two-dimensional stiffness coeffi- 
cients cannot be satisfied by an ordinary shell, in 
which the bending moments result merely from the or- 
dinary stresses distributed in the thickness t, because, 
the bending stiffness coefficients being then proportional 
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to t? and the membrane coefficients simply to t, the 
shell becomes more and more a membrane when the 
thickness tends to zero. But it can be satisfied if the 
shell is a generalized one, |. e., if non-vanishing couple- 
-Stresses are admitted to exist (see [7]). 

Convergence analyses can be based in any case 
on an approximation theorem which states that the dis- 
tance between an arbitrary element s in X, which mini- 
mizes q on €, and its approximation in X', s', , satisfies 
the inequality 


dis <A) + NV Ds pl+Idaç! 


& 


(4.1) 
nhera 

dep=qls)—q(s) (4.2) 
Sap=els,)—gqls,) (4.3) 

and S, is isoconstrained with s and such that 
Sa =B (s;) (4.4) 

Indeed, as s minimizes y on C, 

pis) <gls,) (4.5) 


On the other hand, as s', is the aproximation of 
sin X, s', minimizes q on C”, so that 


pis,)<qls) (4.6) 


Introducing (4.2) and (4.3) into (4.6), there results 


pisglt,p— 8,7 Spls) (4.7) 
Combination of (4.7) and (4.5) vields 
pisgltôp-—-d,p<pls)<pl(s,) (4.8) 
Therefore, 
plsal—pls)<IiS pl+rIid, gl (4.9) 
On the other hand, (2.8) leads to 
dº (ss, )=q(s,)—q(s) (4.10) 
so that inequality (4.9) can be transformed into 
dissa)<NW/ Issgl+Id gl (4.11) 


or, by virtue of the triangular inequality, into (4.1). 
Let us apply the theorem to the case of Ritz's 


method in which C' €C JC. This means that s' and Sa 
belongs 
(4.1) becomes 


ara isoconstrained with s and Sa. Then, as Sa 
to €, Sa 
therefore 


can be identified with Sa 


(4.12) 


On the other hand, as s is isocontrained with s'; 
(2.8) leads to 


dº(ss)=qpl(s)-qis)=ô,g (4.13) 
so that (4.12) vields 
d(s,s' al < dis,s') (4.14) 


The app:oximation error is thus bounded by the 
interpolation e roz. In other words, as far as Ritzs me- 
thod is concerned, completeness implies convergence. 

Inequality (4.1) can be simplified if is is considered 


that T+ is the sum of two terms, one the order 
stiain e 

of the variation modulus - and the 
stress 

| displacement Era 
other of the order of the k | variation 
traction 
ou | : 

modulus RO DR 

ôPp 


sT=0O(Isel)+O(Isul) 8T*=0 (1851) + 


+ O (18pl) (4.15) 


On the othe: hand, by virtue of (2.5), 


diss,)=0O(18, e!) a ss q) =0 (8,51) 
(4.16) 


so that the fist term in (4.1) can be neglected and 
(4.1) be transformed into 


dels s,)<O(ls ul)+oO (la el)+ 


+ O (ls ul) +O(l5,e1) (4.17) 


Any convergence or accuracy analysis [7,8] is thus 
reduced, according to the present theory, to two essential 
steps: a) the construction of a field Sa isoconst.ained 
with s and interpolated by the aproximate solution, 
b) the evaluation of the order of 58qw, i.e. of the 
variation of the functional associated with the interpo- 
lation error, both for the exact and for the aproximate 
solutions. Às this evaluation is not always easy to make 
in connexion with the approximate solution, resorting 
to especial techniques like the so-callad patch-test, be- 
comes sometimes necessary [8]. 


5. Conclusions 


lt would not be fair finishing this paper on the 
applications of Functional Analysis to the Theory of 
Structures without referring Prager [9] and Synge's 
[10] pioneer work, as well as Mikhlin' book [6]. A 


10 


more recent book by Oden [11] should also not ba 
omitted. 

In all such papers no attempt was made however 
to presenting the Theory of Structures, as a scheme of 
Mathematical Physics general enough to become a 
theoretical support not only to discretization but to the 
general problem of the generation of models from other 
models, 

Ths theory of convergence contained in the scheme 
presented by the author can itself still be used outside 
structural analysis, in every case in which a minimum 
principle is known to exist. (see [12]). 
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Programação Linear com Relevo para a Programação 


Linear Inteira e as Variáveis Booleanas 


RESUMO 


Este trabalho, resultante da reunião de dois temas 
por nós tratados no último ano da Licenciatura em Ciên- 
cias Matemáticas, é dedicado acima de tudo à Progra- 
mação Linear. Na sua Parte 1 cinge-se principalmente 
ao estudo, um tanto pormenorizado embora, do Método 
do Simplex e do Método do Simplex Dual tentando 
apoiá-los numa base matemática sólida e simultanea- 
mente aplicá-los a alguns problemas simples de modo 
a completar a sua compreensão. Na Parte 2, aproveitando 
o estudo anterior, dedicamo-nos à Programação Linear 
Inteira debruçando-nos em particular sobre os algoritmos 
da Gomory e de Balas por nos parecerem de grandes 
possibilidades pela enorme gama de problemas que 
podem enfrentar. E a natureza do segundo algoritmo 
levou-nos a um tratamento especial de transformação 
de variáveis inteiras em variáveis booleanas. Programa- 
mos ainda estes algoritmos na linguagem A.P.L. apre- 
sentando de seguida as listagens de alguns problemas 
simples resolvidos pelo computador. 


PREFÁCIO 


O trabalho que se apresenta resulta da reunião de 
dois outros trabalhos de títulos o:iginais «Programação 
Linear com Relevância para o Método do Simplex» e 
«Programação Linear Inteira e as Variáveis Booleanash, 
tratados respectivamente no Seminário e Estágio do 5.º 
ano da licenciatura em Ciências Matemáticas na Univer- 
sidade de Luanda em 74/75. E ocorreu-nos tal reunião 
pelo facto do primeiro ter sido desenvolvido de maneira 
a possibilitar um mais fácil acesso aos tópicos escolhi- 
dos em Programação Linear Inteira. De qualquer modo, 
talvez para respeitar os títulos originais, dividimos este 
trabalho em duas partes, 

Na Parte 1, que engloba os primeiros quatro capi- 
tulos procurámos fazer um estudo tanto quanto possível 
exaustivo do conhecido Método do Simplex e do Méto- 
do do Simplex Dual pela capacidade que estes mé- 
*odos encerram de resolver numerosíssimos problemas 
de aplicação aos mais variados campos e ainda pela 
necessidade de dominar pe:feitamente estes dois méto- 
dos para poder abordar com êxito a Programação Linear 
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JOSÉ DE MATOS LOPES TEIXEIRA 
Matemático 
Assistente do |, S. T. 


ABSTRACT 


This work is a resultant of two former papers 
which have been prepared as a part of the requirements 
for the graduation in mathematics and deals with Linear 
Programming. Part 1 is mainly a study, with certain 
details, of the simplex and the dual simplex methods, 
seen from the point of view of a rigorous mathematical 
approach aiming at the same time to apply them to 
simple problems in order to get a good understanding. 
On the basis of the theory previously developed, Part 2 
deals with Integer Linear Programming, with a special 
enphasis on the algorithms of Gomory and Balas which 
proved to be of particular interest because of their ex- 
treme versatility to attach concrete problems. The nature 
of second algorithm has led to the usage of boolean 
variables instead of integer variables. Both algorithms 
have been coded in A.P.L.; as a part of the work we 
also show listings of the outputs of some easy problems 
solved on a computer, 


Inteira. Foi nosso objectivo basear os algoritmos dos 
Métodos referidos em verdades matemáticas perfeita- 
mente comprovadas e por outro lado apresentá-los es- 
quematicamente no sentido de tornar possível enfrentar 
de imediato alguns problemas cuja resolução apresen- 
tamos com certa riqueza de po:menor. 

Na Parte 2, constituída pelos capítulos 5, 6 e 7, 
enveredou-se então pela Programação Linea: Inteira que 
de resto se desenvolveu à custa das limitações encon- 
tradas no Método do Simplex criando variados algo- 
ritmos. Pareceu-nos contudo aconselhável tratar exclusi- 
vamente os algoritmos de Gomory e de Balas que pro- 
gramámos em linguagem A.P.L. Apresentamos também 
a listagem pelo computador dos sucessivos cálculos 
para cada exemplo tratado e incluímos também um certo 
tratamento da transformação de variáveis inteiras em 
variáveis booleanas. 

Uma palavra de agradecimento sincero para o 
senho: Prof. Dr. J. Marques Henriques que nos orientou 
e entusiasmou a tratar da Programação Linear, Geral e 
Inteira, destacando aqui a importância das variáveis 


booleanas. Lisboa, 14 de Fevereiro de 1975 
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PARTE I 
Capítulo 1 
INTRODUÇÃO A PROGRAMAÇÃO LINEAR 


1.1. — Nascimento e Limites Actuais da Programação 
Linear e não-Linear. 


Quando um responsável por qualquer actividade 
coordena todas as suas tarefas, limitações e interliga- 
ções possíveis da modo a obter o máximo proveito, está 
na realidade perante um problema da Programação Ma- 
temática. Mas para ele ser resolvido com êxito torna-se 
necessário traduzi-lo por expressões matemáticas e de- 
terminar um processo de encontrar o objectivo em vista, 
isto é, formar o modelo matemático que sintetize o pro- 
blema em causa. 

Foram sobretudo as questões de indole econômica 
e militar pela grande quantidade de factores em jogo, 
que fizeram avançar este importante ramo da Matemá- 
tica. Embora já do século XVIII sejam conhecidos mo- 
delos de Programação Linear como o «TABLEAU ECO- 
NOMIQUEs» de Quesnay que tentava equacionar as rela- 
ções de produção entre os propriatários agricolas, foi em 
1936 que Leontief com o seu método da análise das 
tabelas de «input-outputr baseado na ideia de que uma 
grande proporção do esforço da Economia Moderna se 
dedica à produção de produtos intermédios e que a 
produção destes produtos intermédios tam por sua vez 
grande influência na produção de produtos finais, que 
se deu o verdadeiro arranque. 


Mas o criador da Programação Linear foi George 
B. Dantzig que em 1947 apresentou uma técnica para 
planificar as diversas actividades da Força Aérea dos 
Estados Unidos. Esta técnica tem bastantes semelhan- 
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PARTE II 


Capitulo 5 Programação Linear Inteira. 
5.1. Introdução. 


5.2. Aplicações e Modelos. 


Algoritmo de Gomory e sua Programa- 
ção em A.PLL, 


Capítulo 6 


6.1. Generalidades. 

6.2. Escolha do Plano de Corte. 

6.3 Selecção de A. 

6.4.  Esquematização do Algoritmo. 

6.5. Exemplo Resolvido. 

6.6. Programação do Algoritmo da Gomoy. 


Variáveis Booleanas e Algoritmo de 


Balas. 


Capitulo 7 


7.14. Transformação de Variáveis Inteiras em 
Variáveis Booleanas. 

7.2. Métodos de Enumeração Implicita. 

7.3. Algoritmo de Balas. 

7.4. Exemplo Resolvido. 

7.5. Programação do Algoritmo de Balas em 
APL, 


ças com o método de Leontisf pois analisa as relações 
que as dive:sas actividades de rectaguarda têm com o 
objactivo final, relações análogas às existentes entre os 
produtos finais e a produção dos sectores intermédios 
no método de Leontief. Mas a inovação de Dantzig foi 
consegui: um método para decidir qual dos vários 
planos que satisfazem os objectivos finais é o melhor, 
portanto um método de optimização. Também em 1947, 
T. Koopmans desenvolveu estes métodos aplicando-os 
ao chamado «Problema de Transportes. 

A Programação Linear, como o nome indica, debru- 
ca-sa sobre problemas em que as suas restrições ou 
relações são da forma: a, x, tax, +... ta, x, =b, 
onde os a, e b são coeficientes conhecidos e os x, são 
variáveis a calcular, 

Parece que foi Tourier o primeiro a estudar siste- 
mas de desigualdades lineares tendo realçado a impor- 
tância desse estudo na Teoria das Probalidades. 

O verdadeiro início do cálculo da Programação 
Linear comecu com a publicação do «Método do Sim- 
plox» por Dantzig de parceria com Leonid Hurwicz e 
Koopmans. Actualmente as aplicações de Programação 
Linear estendem-se a todas as categorias de actividades 
industriais, comerciais e governamentais e a versatili- 
dade e capacidade do modelo de Programação Linear 
parece não ter limites. Mas até os primeiros investiga- 
dores lhe reconheceram certas limitações e assim se 
começaram a desenvolver bases teóricas e processos de 
cálculo para certos problemas especiais de Programação 
não-Linear, isto é, problemas onde as restrições já não 
são do tipo anterior, rectas, planos ou hiper-planos, co- 
mo por exemplo: 
minimize a função f(x, x,) = — log x, — log x, sujeita 
à restrição x, + x, < 2, onde x, e x, são logicamente 
não-negativos. 
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A partir de 1950 um grande número de investiga- 
dores (matemáticos e economistas) dedicaram-se ao 
desenvolvimento das diferentes ramificações da Progra- 
mação Linear: Programação Quadrática, Programação 
Estocástica, Programação Separável e Programação 
Linea: com variáveis inteiras onde desempenham um 
grande papel as variáveis booleanas nos métodos de 
«Enumeração Implicita»s. Os domínios de aplicação 
multiplicam-se e a dimensão dos problemas resolúveis 
aumenta rapidamente provocando o interesse crescente 
dos especialistas de organismos públicos e privados e 
das Empresas Industriais e Comerciais, Claro que muito 
deste interesse provém da utilização de computado-es 
cada vez mais pode-osos permitindo a resolução de 
problemas de até vários milhares de restrições. 

As possibilidades abertas pelos problemas de Pro- 
gramação Linear com variáveis inteiras (P.L.l.) ou «mis- 
tas» (isto é, problemas onde todas ou uma parte das 
variáveis devem ser obrigatoriamente inteiras), são de 
tal modo vastas que várias equipas americanas e euro- 
peias a eles se dedicam na preparação e desenvolvi- 
mento de algoritmos cada vez mais eficientes. 


1.2. — Considerações gerais e Notações. 


Num problema de Programação Linea; procura-se 
uma repartição mais eficaz dos recursos disponíveis 
tendo em vista certo objectivo. Imaginemos uma em- 
presa que se dedica às actividades X,, Xy, ..., X, + 85- 
tando os investimentos, as despesas de armazenagem, 
as horas de trabalho, o consumo de energia limitados 
respectivamente às quantidades ou valores b,, b,, b,, 
Db, -.. etc. Podem determinar-se na empresa os in- 
vestimentos 8,;+ 05 custos de armazenagem a,,, os 


tempos gastos a,;, , os consumos de energia Ag ques 
etc., para cada unidade produzida pelas diversas activi- 
dades x, . Depois de calculados os lucros unitários c; 
para cada actividade, podemos então formular o pro- 


blema: 


Maximizar o lucro Z=€CX TCX Te. TC x, 


sujeito às restrições 
dai Bem Teco Ta, A <b 
Ao, X; T Boy Xy To TB, X <b. 


—— ss = (= E == 


* - = 
Da miXk TAmak Too tan Xp Eb 


e COM Xjs Xp «uu XX, >0 


Suponha-ss agora que um ourives produz dois tipos 
de anéis, A e B, usando ouro e prata em doses variá- 
veis. Consome três gramas de prata e um de ouro no 
tipo A; dois de prata e dois de ouro no tipo B. Sabendo 
que só possui 52 grs. de prata e 28 grs. de ouro e 
que além disso o lucro por anel do tipo A é de três 
unidades e do tipo B é de 4 unidades, calcular qual 
o número ideal de cada tipo de anéis que o ourives 
deve fazer de modo a maximizar o seu lucro. 

Sejam x, & x, Os números de anéis a produzir, res- 
pectivamente dos tipos À e B. Teremos 
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Maximize L=3x, + 4x, 


sujeita a 


x + 2x,4 28 6 xx, 40 


3x, + 2x, < 52 


cuja solução, determinada pelo Método do Simplex (ver 
4.2., 6.º) 6 x,=12 e x,=8 com Z=3X 12+4X 
xXx 8B= 98 


Podemos então dizer que em Programação Linear 
o problema consiste na procura do vector X = (x,, Xx,, 
+ X, | que extreme uma função 


UX) =c,x +6% +... toe x, 
de valor real definida no espaço n-dimensional e linear, 
a qual recebe o nome de função objectivo, e que simul- 
taneamente satisfaça as restrições igualmente lienares: 


X; > O para j=1,2,..an (1.1) 
e 21X FagpX to ta, xp =b, 
Bo, X, TA XT a, Xp, =D, 
do di a a (1.2) 
mi Tank To tax =D, 
em que o vector linha c=(c,C sc, )JER". 


cujas componentes são chamadas coeficientes de custo, 
o vector coluna b = (b,,b,,...b JER” 
nentes a, são constantes reais dadas (com m < n) 


8 às compo- 


Po: vezes escreve:-se-á apenas: 
[8] n 
PR as + sã 
Minimize > Cc; x; sujeita a 3) apx; =b;, x,>0 
Il= me 


s Comit= 12 as mboli=A 2 cn) 
ou ainda mais sinteticamente 
Minimize cX sujeita a AX=b e X50 


onde X e O são vectores coluna n-dimensionais, O com 
componentes todas nulas. 


NOTA: Para facilidade de expressão substituir-se-ão 
por vezes as componentes a, da matriz A por 
X;; + OS vectores coluna da mesma matriz por 
P. eo vector coluna b por P, , isto é, a res- 
trição AX = b passará a ser substituída por: 
MP trxPt.. Tx PSP, 
Convém agora recordar algumas definições e teo- 
remas de grande importância na sequência deste tra- 
balho: 
Definição 1 Chama-se solução admissível dum pro- 
blema de Programação Linear ao vector 


X = (X Xor ...,X, ) que satisfaça as con- 
dições (1.1) e (1.2). 
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Definição 2 


Definição 3 


Definição 4 


Definição 5 


Detinição 6 


Definição 7 


Definição 8 


Definição 9 


Teorema 1 


Teorema 2 


Teorema 3 


Teorema 4 
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Chama-se solução básica admissivel à 
solução de componentes não negativas 
obtida, fazendo n-m variáveis iguais a 
zero, para as m variáveis restantes, quan- 
do os seus coeficientes não formem de- 
terminante nulo. 


Chama-se solução optimal à solução ad- 
missiível que minimize ou maximize a fun- 
ção objectivo. 


Chama-se combinação convexa dos pon- 
tos Yy Yo Y,y E EM ao ponto Y - 

q Yi tayat... +a,y, coma, es- 
calares não negativos e tais que 


Um conjunto CE EM é convexo se e só se 
para todos os pares de pontos y, ey, a 
sua combinação convexa y=a,y, + 
+ a, Y, também pertence a esse conjunto. 


Um ponto dum conjunto convexo diz-se 
ponto extremo se não for possivel expri- 
mi-lo como combinação convexa de quais- 
quer outros dois pontos distintos. 


Chama-se ambiente convexo C(S) da 
qualquer conjunto S dado de pontos, ao 
conjunto de todas as combinações con- 
vexas dos pontos de 5. 


Chama-se poliedro convexo ao ambiente 
convexo do conjunto S formado por um 
número finito de pontos. 


Um Simplex é um poliedro no espaço a 
n-dimensões e que tem exactamente n + 1 
vértices. Às fronteiras de um Simplex 
contêm ainda outros Simplex de menor 
dimensão 


O conjunto das soluções admissíveis dum 
problema de Programação Linear é fecha- 
do, éu M=(X:XER,X5>0,AX=5) 
é fechado. 


O conjunto das soluções admissíveis é 
um conjunto convexo e tem um número 
finito de pontos extremos. 


A função objectivo alcanca o seu minimo 
num canto do conjunto convexo M forma- 
do por todas as soluções admissíveis. 


Se M é não vazio e limitado, a função 
objectivo tem minimo e assume-o em pelo 
menos um dos seus cantos. 


Teorema 5 Se M não é limitado mas se a função ob- 
jectivo possui minimo, esse mínimo é 
atingido pelo menos num dos seus cantos. 
Teorema 6 A solução X = (x, x, ...,X, ) é um can- 
to de M se e só se as x, componentes 
positivas são coeficientes de vectores li- 
nesrmente independentes P, em número 
de m em: 

dk Pt Pd PP =P 


n 


Ls 


NOTA: Estes teoremas são de importância fundamental; 
contudo a fim de não alongar demasiado este 
trabalho não são apresentadas aqui as demons- 
trações que podem ser lidas no livro de Saul J. 
Gass, Programacion Lineal, Capítulo 3. 


Analisando bem o conteúdo destas definições e 
teoremas podemos seguir um raciocínio semelhante 
aquela que levou G. B. Dantzig ao Método do Simplex 
a que vai ser motivo de particular interesse neste tra- 
balho: 


a) A função objectivo toma o seu minimo num 
pelo menos dos cantos de M. 


b) A cada canto de M estão associados m vectores 
lingarmente independentes do conjunto dos n vecto-es 
dados. 


c) As soluções do problema são cantos e há tan- 


tos como ( b- | maneiras de obter conjuntos de m vec- 
tores linearmente independentes. 


d) Sem e n são paquenos não haverá muitas so- 
luções a analisar mas o vulgar é serem m e n números 
bastante grandes o que obriga a seguir um método que 
da algum modo nos leve à solução óptima. Um desses 
métodos e o mais eficaz é o Método do Simplex que 
consiste sumariamente no seguinta: 


1.º Determina-se uma solução de ponto extremo a 
partir de um conjunto de m vectores linearmente inde- 
pandentes e, por aplicação do teorema B, (ver 2.2.2), 
verifica-se se esta solução minimiza a função objectivo. 


2.º Se não a minimiza determina-se uma outra so- 
lução de ponto extremo vizinho, isto é, ligado ao ante- 
rior por uma fronteira do poliegdro convexo, e que dê 
à função objectivo um valor inferior ao anterior. 


3.º O trabalho da alinea anterior repete-se um certo 
número de vezes e geralmente entre m e 2m vezes. 
Termina o processo quando se atinge o mínimo. 


Capítulo 2 
DESENVOLVIMENTO DO METODO DO SIMPLEX 
2.1 — Solução Admissivel Inicial. 
Ao ser-nos dado um problema de Programação 


Linear as suas restrições podem apresentar um dos 
quatro seguintes aspectos: 
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1.º O modelo matemático é da forma 
AX=b com b50 e m<n 


em que se supõe que pelo menos uma b, >0 e temos 
as duas situações seguintes: 


a) A matriz À não contém vectores colunas unitá- 
rias e então inicia-se o cálculo com m vectores artifi- 
ciais, técnica esta chamada de «base artificial». Ora o 
problema geral de Programação linear é: 


minimiza CM TOS Te FO x. 

sujeita a anXt..ta,x=b 6 x>0 
dnX To Tax =b, 
da T+ Anta Sb, 


e então dentro desta técnica aumentamos o sistema 
pa:a: 


minimize 
GX% +. + Xy EWXn4 FWXnyo + WX n+m 
sujeita a 
an X Te Tan XT X+1 mio 
Bite E es E ga Mi TF Xn+9 =b, 
das Es FB Ha E Xnymo Dm 


e ainda X; >0 para j=1,2, ..,n+m. 


em que w representa qualquer quantidade suficiente- 
mente grande que nada altera por as variáveis corres- 
pondentes acabarem por se anularem. Assim os m vec- 
tores Pig Poço cs P,sm formam uma base (arti- 
ficial) para este problema aumentado e se o problema 
inicial tem solução óptima, o problema ampliado tam- 
bém a terá e a própria estrutura do Método do Simplex 
garante que seja impossível aparecer uma das variáveis 
X,4j com um valor positivo. Se o problema original 
não é possível, então a solução do problema aumentado 
apresenta:á uma variável Xn4i 20. 


b) A matriz À contém k vectores unitários distintos 
a então acrescentar-se-ão m — k vectores artificiais e 
aplica-se o número anterior. 


2.º O modelo matemático é da forma 


AX <b com b>O0 

Neste caso transformamos cada desigualdade numa 
igualdade acrescentando a cada linha um variável não 
negativa chemada «variável da folgar e que faz apa:ecer 
uma base de m sectores unitários P,44+ Pior 
P para iniciar o processo: 


e Ê 


n+m' 


anXM rt... tax &b, 


81X +... + agnxn SD: 


Em X, Fe + êmn kn Sb 
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passa a 
2%, Ft Aga + Xna4 = b, 
a, X É... Tt A9nXn E Raso =b, 
Dik T+: TT Amado E Xntn = D 


3.º Se fo: da forma 
AX > O com b>0 


podemos igualmente transformar cada inequação em 
igualdade como na alínea anterior mas desta vez apare- 
ce-nos um conjunto de m vectores unitários negativos, 
ou seja AX—X=b com X=(XgeXnto cce Xnsm) 
um vector coluna não negativo, o que nos obriga a um 
segundo artifício para encontrarmos os m vectores uni- 
tários positivos. Transformamos AX-X=b da se- 
guinto maneira: 


a; =-—at+aç--- para [ESLjOKZ e rm 
a ã,j j=1,2,...n+m 
bi = —b +b...parai “se sendos o índice de max b, 
b';=b, 


ou seja, determinamos max b, =b, e somamos a 
s-ésima fila à simétrica de cada uma das outras o que 
fará aparecer apenas m — 1 vectores unitários positi- 
vos e que serão completados com um vector artificial 
correspondente à s-ésima variável pelo processo 1.º a). 


Vejamos o exemplo: 


* =X Pk al Meio = 1 
MPE MT Ah” =k = 0— 
dd Fte — Ri 254 + Mm dd ca =x = 

2x, — 4x, + X, =%4=1 
= =x 2 — Em +Xx—x=2 
MO ta —X = 2 


a aplicávamos ao último sistema o processo de 1.º a). 


4.º Um último caso possivel é aquelas em que apa- 
recem misturadas equações e tipos diferentes de desi- 
gualdades: 

AX = b com b sem qualquer restrição 


onde por manipulações apropriadas dos números ante- 
riores se podem fazer surgir os m vectores unitários. 


2.2. — Geração de Soluções de Ponto Extremo e Pro- 
cura da Solução Minima. 


2.2.1. — Geração de Soluções. 


Após a aplicação de algumas das alíneas ante- 
riores, O que não é necessário se a matriz do problema 
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apresentar em evidência o conjunto de m vectores linear- 
mente independentes podemos iniciar o método: 

Suponhamos que o conjunto dos m vectores P, 
linearmente independentes são os primeiios m. Então 
o vector solução inicial, ponto extremo, será: 


X = (XX co X + O, ..., O) com x; O e 


m 


Bt ERP co Tx Pr=P 


HH NR (O 


(2.1) 


como nos garantem os teoremas anteriormente enuncia- 
dos. Sendo P,, ..., ma linearmente independentes pode- 
mos exprimir como combinação linear deles, qualquer 
um dos n vectores dados. 


m 
Bo = a N ; = 4 
P, Pã Xj P, para j= 1,2,.. un 
j==1 
Be onde Xj =aj conforme ficou assente na nota de 


(1.2). Suponhamos agora que P. com k =m é um dos 
vectores que não se encontra na base e que tem pelo 
menos uma componente x, >0 na expressão: 


P (2.2) 


REG + AP Fo Ra Pa 


Multiplicando esta igualdade por um número qual- 
quer 9 e subtraindo-a membro a membro de 2.1), vem: 


e; — EMP, + iã— OnI)Ps, Fo + 
+ (Xm 0X mk | Pt Pg =P 
e o vector X =(X, — 6X, pr Xo — 8X porra Xy— 
E Mn 


é uma solução possível se todas as suas com- 
ponentes fo-em não negativas, o que sucede se fizermos 
8>0, que ainda obriga a ser X'>X. Como algumas das 
Xik podem ser não positivas vamos seleccionar 4 > O 
tal que 


XxX. —8X Ike Ge 


O ou 0<x; [X; para todas as Xi >0 


Ora qualque: 9 que satisfaça 


x: 
o 4 
O<e<min — 


I Kik 


da:ia uma nova solução possível X' com m + 1 com- 
ponentes positivas. Mas como queremos atingir uma 
nova solução de ponto extremo, atendendo ao teorema 6, 
não podemos ter mais de m componentes positivas, e 
teremos de escolher um 6 dentro daquele inte:valo mas 
que anule uma das m + 1 componentes. E se escolher- 
mos precisamente 
E ei Xj 
É — min RR 024 =— 
i x k 

1ê-seo imediatamente que se anula a componente de X' 
fornecedora desse mínimo. Admitamos que ela é a pri- 


meira componente: 6, = x,/x, + então fazendo 
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— 69 Xi para = ssm) O 


n 
k 


1] 


Xm+1-%, além de Pr= Psi 


temos 


X Pa + xs dé RU Xm Pm + X m+1 Pm+1 ci Ps 

Claio que se todas as Xi k fossem não positivas, 
não conseguiamos determinar um 4 que anulasse uma 
das m + 1 componentes e obteriamos apenas uma so- 
lução possível associada com os vectores P,, P,, ..., 
P,n41 8 isto indicaria que o problema não teria solução 
mínima finita como veremos em estudos seguintes. Mas 
para se afirmar que X'= (XX. X'm+4] é ponto 
extremo ainda falta provar: «os vectores associados 
Por Par coor Pm+.q São linearmente independentes»: 


Demonstração: Admitamos que eram linearmente 
dependentes. Então 


dy Po + dy Pa ee + dm Pm + dm Pm+1 0 


onde nem todos os números d, são nulos. Ora P,, P., 

- P São linearmente independentes po:que é um sub- 
conjunto de P,, P,, ..., P,, linearmente independentes 
por hipótese e isto obriga a ser d,.,., “ O. Assim, pondo 


E = — d [dn+4- virá a igualdade: 


Ped 6 PF TP aa BP 


m m om 


que subtraída da P 44 =P =Xk P +- X k Pa unisk 


+ x P 


“ mk k 


anterior dá: 


X,kP; = (x, k — 8.) Fo H- (Xa k nc e.) E + 


o a É = Bad ra 4 


m,k 
em que todos os coeficientes dos vectores terão de ser 
nulos por serem estes independentes, o que não sucede 
para x, que se supõe positiva. Logo concluímos que P,, 
Pas... P, São linearmente independentes e X' é solução 
de ponto extremo, 

Assim com o vector P, se encontra em função 
da nova base, temos de exprimir também todos os ou- 
tros vectores P,, P.,..., P, em função da nova base 
que já não inclui o vector P,. Ora de (2.2), ou saja 


Ei == Mk P, à Ka k E "Te =. + Xm.k Fa” sai 
P, se P E Ko k ”. ns ia, P ] 
Xi, k à k m,k * m 


Então para qualquer vector P. que não esteja na base 


original, ou seja 
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a 
mi a E Kay) Pa 
Tik 


a 


E RR É 


paia (j;-k, 2, 3, ..., m), semelhante a P, =6P, + 
+ (x, o 0X. k) Po + (x, — 0X ) P, Paes 


Agora dispomos duma outra solução de ponto ex- 
tremo associada a outra base e em função da qual estão 
expressos todos os outros vectores, Podemos novamente 
fazer entrar na base um outro vector e obter outra base 
e nova solução e assim sucessivamente até à solução 
óptima se ela existir. 

A importância dos coeficientes da expressão de P 
anterior justifica a seguinte: 


NOTA: Vê-se que é passagem fundamental do Método 
do Simplex ir exprimindo todos os vectores P, 
em função da nova base. Suponhamos então que 
a base inicial é (P,, Pa... Pal=1, Ivec 
tores unitários). Podemos exprimir todos os ou- 
tros vectores em termos desta base, que serão: 


Po = XoP, Ts + X9P| ce Ma K” 


P — Rd E. + > UE + XP) = ea é + Kmk Eli 


i m 


P, > Mi E: Eos FX P, Meios PF Xm. F 


Mas se agora se concluir que deve entrar na base o 


vector P, e sair ovector P, , da expressão de P,., 
temos: 
1 X X 
1,k ak 
P, = pp —..— É po 
x A | 
Xi x |,k L.k 
e substituindo este vector P, em P, .... P,. gotas) P, 


- P, - vectores não incluídos na nova base, vem: 


*o Po | E 
Po “(Mo Mask) P, + ixo —— Xa k)Po a 


1,k *Lo 
iza X 
lo lo 
| - P,. +... + mo — co Xmk) Em 
Xk Lk 


(continua na coluna seguinte ) 


n.º da base  Coefici. Solução. 
linha inicial básicos 


inicial 

B C Ps 
1 P, Cc, *1,0 

P, C. As 0 
3 P, Cs AZ o 
Em Em Em Xm,o 
m + 1 | Eá 
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» Mi “1 
Pa EO By PST OG — = Sal; 
*1k XI k 
desci É sao di À as P 
j E mi Xmk) Pm 
uk E 


em que P, é a nova solução (vector coluna) e P, os 
outros vectores coluna do novo quadro (veja-se exem- 
plo final), E temos então; 


X kj > 
Xk 
que são chamadas «FÓRMULAS DE ELIMINAÇÃO COM- 
PLETA» que nos permitem obter a nova solução e as 
novas representações dos vectores que não se encon- 
tram na base. 


2.2.2 — Procura da Solução Minima. 


Temos então uma solução inicial do Mio co 
Lis “4 
o Xmol & O correspondente conjunto de vectores 
linearmente independentes P,, P,,...,. PP, isto é 


Kofi tMoPa Tt +rmoPm=P 


mo m É] (2.3) 


* Paio t *2,0 Co Te E Xm,o Co — Z, ou Poe =Z, 
(2.4) 


onde Z, é o valor da função objectivo, C,, C,, ..., Cy 
são os coeficientes de custo e X, 920. Expressando cada 
vector P. em termos da base P,, P,,.... P, OU 


P=x5P,+xaPa +... + xmyPmjcomij=1,2,...n, 
[2.5] 
definimos então 


2 = M10 + %]96 To Tx 


e. mi= dA ash 


mo “m 
(2.6) 


onde cada c, é 0 coeficiente de cada P, da base. 

Para melhor compreensão de notações e melhor 
transposição para a prática, apresentada no final do 
trabalho, convém observar o seguinte quadro (seme- 
lhante aos usados na prática manual): 


zeros ou w (coefici- 
entes artficiais) 


Pn+1 Pntm+1 

1 O O 

[8] O (8) 

[6] 1 O 
tou (—1)) 

(6) OQ 1 

O Q O 
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e assim podemos enfrentar melhor os seguintes teore- 
mas: 


Teorema 7: Se Z; —c >0 para alguma coluna j, 
então existem outras soluções tais que os 
correspondentes valores Z da função objec- 
tiva são tais que Z<Z, , onde o limite 
inferior de Z pode ser finito ou infinito. 


Demonstração: Como vimos, multiplicando a igual- 
dade (2.5) po: 4 e subtraindo-a de (2.3) obtemos uma 
nova solução de ponto extremo no caso de alguma 
X;j >0 em (2.5), pois assim 6 passa a 6, =min (Xjo/ 

| 
[Xi | que anula uma das m + 1 componentes da nova 
solução 
Pi= ix “6X, |) P 


1,0 E eco Pag Ma VE PP 


cujo correspondente valor Z da função objectivo será: 


Mo Mg) + Fo” Xmj! Cm TOC = 
=*%,0€% To TXmo Cm” 
— 604; CT... + Xl Toc = 
= 8 Tt =, IZ —E |) 
ou seja 
£=Ey-U(Z, —e;) 
mas como 9>0 e por hipótese £ —€ >0, virá 


2<£, - Claro que nas condições postas, para alguma 
X,,20, 6 = 6, será finito e portanto Z terá um limite 
inferior finito; se pelo contrário todas as Xi <O, então 6 
não terá limite superior e portanto Z > — O e neste 
caso para qualquer 9>0 obtém-se uma solução com 
m + 1 componentes e logo não extrema. 

Mas o que nos interessa é que sendo Xi do vec- 
tor coluna P, a entar na base, maior que zero, se ob- 
tém nova solução de ponto extremo e nova base, à qual 
corresponde Edo - Esta nova base pode ser usada 
como a anterior e se novamente existir um valor 
Z; — e; >0 e uma cor:espondente Xi >0, encontrar- 
-se-á para a função objectivo um valor Z'<Z e assim 
sucessivamente. 

Quando termina o processo? O teorema seguinte 
responde a esta pergunta e isso verifica-se quando to- 
dos os valores da m + 1 — ésima linha, Z —c, |, fo- 
rem não positivos ou quando para algum Z, ade 
se tenha X;j < oO. 

Teorema 8 Quando for £ —c <O para j=1,2, 
“s n, então a solução de ponto extremo 
X=(Xo cc Xmo) é também solução 
minima, isto é, PG =XoP,+tXo Ps,+ 
TF ok Po 02 o GT tio CG 


Fo + Xmofm= Minimo. 


Demostração: Vamos tomar uma outra solução possivel 
para a qual o valor da função objectivo seja Z' e verifi- 
car que Z, < Z', Seja então a solução 


Po=Yo PiTYrag PT TYho 
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2 = y,08 + Yagla Te. Ty 


ER C, 


Uma vez que Z — €; SO para Jj=1,... h OU 
Z; <e; . substituindo c, por Z. na igualdade ante- 
rior, vem: 


L >Y,0%: ks Y. 0fs +Y Z, 


RE! 


Substituindo nesta desigualdade, a igualdade de defini- 
ção de Z, ou seja 


e reordenando os termos pelas propriedades comutativa 
e associativa, vém: 


n mn 
+ É a, ” ' 
2> (E %on)a + (Bo xi jato + 


951 i= 


in 
ApY% y, 
| -— Tio Km, Em 


| 


Basta prova; agora que o segundo membro anterior é Z,, 
Ora pela definição desta quantidade, é preciso que os 
conteúdos dos parêntesis sejam simultaneamente os 


coeficientes da nova solução P, e iguais a X, 9. X 9: 


i Xm.o - Para isso façamos a PS di Yo P, + Yao P, Ei 
TE ea ZTE Em ara um tratamento semelhante ao feito a Z.. 
Substituindo cada P. por P =x, P +... + 
sã *mi Em + V8M: 

na o! 

E tm A a Ma 
á “ho doi Mil “ii + tal dad in P, ) nda 

=: = 


Como os vectores P,, P,, ... P, são linearmente in- 
dependentes, os seus coeficientes terão de ser iguais 
aos de P,=XoPitXoPs tr... +Xmo Pm é por 
tanto 


FT Xa00 Te TX Ga É 


Mo cm E] 


ou Z' > Z, O que termina a demonstração. Assim estes 
dois últimos teoremas permitem ir gerando soluções 
admissíveis de ponto extremo e convergindo para a 
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solução óptima (mínima) no caso dela existir ou então 
afirmar que a solução óptima não existe. 


2.3 — Fenómeno de Degenerescência e Técnica de Saida 


No decurso do estudo anterior do Método do Sim- 
plex suposemos sempre a existência de soluções de 
ponto extremo não degenaradas, isto é, dada a base (P,, 
Por cs Po ). não haveria nenhuma componente x, 
nula na solução 


O 


toras Po Tc FMs Pa SP 
e isto para ga;antir que a solução seguinte seja inferior 
à anterior e assim atingir o minimo. De facto se uma 
Xio = O, então viria 6 = ia lxo/x%;1=0 e, pelo 


teo:ema 7, da igualdade Z=Z, —91Z, — G; ), tor- 


nar-se-ia impossível conseguir Z<Z, 

Mas sucede com frequência aparece:-em soluções 
degeneradas e assim ao escolher 6 = O para certa x, ,, 
fazemos entrar na base o vector P, e aparecerá uma 
nova solução básica mas o valor da função objectivo 
permanece constante. Contudo no passo seguinte pode 
já não aparecer nenhuma componente XtsT O e o méto- 
do leva ainda à solução. O problema surge quando há 
mais de uma x, , = O, e portanto dúvida na escolha do 
vector a sair da base, e então aqui podemos cair numa 
situação de periodicidade em que as soluções se repetem 
ao fim dum certo número de passos. O mesmo sucede 
quando embora não haja X; qnulas, se apresente uma 
situação de empate ao obtermos o mesmo 4 para duas 
X io diferentes, pois no passo seguinte aparecerão ne- 
cessa:iamente X,,, nulas. 

Ambas as situações, conquanto pouco frequentes 
na prática são contornadas seguindo-se a chamada 
«Técnica da perturbação» exposta no livro de Saul J. 
Gass e que se resume ao seguinte: 

Suponhamos que o vector a entrar na base é PP... 
que a base é (P,, P,, ..., P, ) e que temos 


indecisão quanto ao vector a sair (P, ou P,?). Então 
calculamos, começando por | = 1, as relações x IJ [Xp 
para todas as linhas da coluna 1 onde houve repetição. 
So ainda apalecerem empates fazemos o mesmo para 
J+1, j+2,... até ao seu desaparecimento. Assim 
neste exemplo calculamos 


At 
min ——= para 
i Xi, 


|= 1,2 


se esse mínimo for x,,,/X,  elimina-se da base o vec- 
tor P,; se for x,,,/x, , elimina-se P,. Em qualquer caso 
é sempre o vector P, que entra na nova base, Se por 
acaso X,1/X,; p =X,,1/X, p + formavam-se as relações 
Xa/Xk O Xya (Xp 8 assim sucessivamente até não 
haver dúvidas quanto ao vector a sair. 
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2.4 — Referências ao Método do Simplex Revisto 


Confo:me já se viu, o cálculo que nos permite ir 
passando de uma solução de ponto extremo para outra, 
é expressar os vectores que não se encontram numa 
certa base em termos dessa mesma base e isto com a 
ajuda das «fórmulas de eliminação completas. Em se- 
guida calculamos os elementos da m + 1 — ésima linha, 
Z, — C, para ver que vector deve entrar na base ou 
ver se a solução já é óptima e determinamos o vector a 
sair da base. Nesta altura repetimos o ciclo entrando 
com as fórmulas de eliminação completa. 

No entanto há um outro método, chamado Método 
do Simplex Revisto que nada acrescenta a este em 
rigor e possibilidades, razões pelas quais não será aqui 
exposto mas que tem as seguintes vantagens: 


a) Em problemas de matriz contendo muitos ele- 
mentos nulos, o número de passagens é reduzido e por- 
tanto o tempo de cálculo é inferior, 


b) Os elementos não nulos da matriz originária 
podem ser armazenados mais compactamente na me- 
mória do computador. 


c) No método do Simplex inicial tem de registar- 
-se uma tábua completa em cada uma das passagens, 
mas usando este, isso não se torna necessário pois basta 
registar apenas uma certa mat:iz inve:sa e O vector so- 
lução. 

Estas vantagens resultam do facto, e aqui reside 
a diferença entre este método e o estudado, de que em 
vez de transformarmos todos os elementos das sucessi- 
vas tábuas do Simplex, transformamos apenas os ele- 
mentos de uma matriz inversa (m X m), sempre dife- 
rente em cada passo, por meio das fórmulas de elimi- 
nação completa. 


Capítulo 3 


DUALIDADE 
3.1. —- Generalidades 


Para cada problema de programação linear da for- 
ma que temos tratado e a que passaremos a chamar 
Primal, existe um outro com ele relacionado e a que 
se chama Dual. E por definição os problemas primais- 
-duais dividem-se em dois tipos: 


1.º — Primal - Dual assimétricos 
Primal —> min f(X) = cX com AX = be X>0 
Dual -—> max g(Y) = Yb com YA « ce Y qualquer 


2.º — Primal - Dual simétricos 
Primal > min fiX) = cX com AX >b e X>0 
Dual -—»max g(Y)=Yb com YA<ceY-o 
em que 
XE TU VEM rca YE 
O, vector coluna n-dimensional 


m! 


A = (a, ) com l=1,2 om: j=1,2...,nem<n 
| 4 Db= ba Dy cos Bo) 


C— (O; er o, É mi 


- n 
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Estes dois tipos podem resumir-se apenas a um 
pois veremos que os problemas primal-dual simétricos 
se podem apresentar na forma assimétrica. Ver-se-á 


também que a variável Y do Dual será Y = c'B”! quan- 


do a solução do Primal é min f(X) = c'X' em que B”! 


é a inversa da matriz (P,, P,, ..., P, ), base da solução 


m 
final. 


Nota-se logo que no Primal o produto AX dá um 


quad;o com m linhas e n colunas enquanto que em YA 


aparecem n linhas e m colunas já que YA — AT yT,o 
que apresenta vantagens pois dada a interligação que 
veremos entre as soluções óptimas dos dois problemas 
pode suceder que o Primal ultrapasse a capacidade do 
computador pelo número de linhas ou colunas e que na 
forma Dual isso se não verifique, ou vice-versa. 


Recordando agora os problemas expostos no início 
do trabalho, é fácil dar uma interpretação económica 
do problema primal e há também uma interpretação se- 
melhante para o Dual se bem que não seja muito evi- 
dente o significado da função objectivo e o das restri- 
ções, contudo, transformando o problema Piimal-Dual 


simétrico para: 
Primal —> max f(X) = cX com AX < be X >0 
Dual —>min g(Y)=Yb com YADceYr>o0 


podemos ap.esentar as formas verbais seguintes: 
Problema Primal: Quantas unidades de cada produto 


(x, | se devem produzir de modo a maximizar o valor 
da produção, sendo dados o valor unitário (c, ) de 
cada produto e um limite superior para a disponibilida- 
de de cada recurso b,. 


Problema Dual: Que custos ou valores unitários se de- 
vem atribuir a cada recurso (3; |) de modo a minimizar 
o valor total dos recursos gastos sendo dados uma dis- 
ponibilidade de cada recurso (b, | e um limite inferior 
para o valor unitário de cada produto tc; |), ou mais 


concretamente: 


Suponhamos que um industrial, por motivos de 
concor:ência, deseja vender mais barato um certo pro- 
duto. Para manter uma certa percentagem de lucro, veri- 
fica que para isso esse produto teria de sair da linha 
de produção a c (valor minimo do custo do produto). 
Então necessita sabe: que quantidades teria de gastar 
dos seus recursos fixos (homens, máquinas e/ou ma- 
térias primas) ou a que preços unitários (93 |) lhe fi- 
cariam esses recursos de modo a gastar o minimo de 
recursos, ou mais simples: Querendo produzir a tanto 
(menos), a quanto lhe vão ficar os recursos utilizados 
(mínimo possível)? No caso de serem homens os re- 
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cursos (por exemplo), qual o seu preço virtual? Claro 
que o industrial, pela diferença, mais valia, entre o or- 
denado pago e o preço virtual vê logo quanto ganha 
ou perde. 


Conforme já se viu pelo estudo anterior, qualquer 
forma de problema de programação linear pode ser alte- 
rado pela junção de variáveis de folga, variáveis artifi- 
ciais e pela igualdade min f = — max (— f). Assim, po- 
dernos transformar os problemas simétricos em proble- 
mas assimétricos: 
simétricos 

Primal —> min f(X) = cX com AX = be X >0 


Dual —Smax g(Y])=Yb com YA< ce Y=>0 


Ora o problema primal transforma-se facilmente no equi- 
valente assimétrico transformando apenas AX -b em 
AX = b com as m variáveis de folga não negativas 


Z = (z,, Z, «-., Z |, vector coluna, pondo: 
(C | O) — C,, l., esp Co D,. D., = nd Dm 
4 dn”! ú 
AI A = 
O au 


Emo à mn 


X Ê, 
[= ) = (Xi Rg a Rs E Ear con E po 


que dá; 


EIA) = (c10)( m ) com (a t-n(s J=» 


B Mx Y 0 
forma primal do problema assimétrico. 


Ora o Dual deste Primal, dual do problema prima! 
assimétrico: max g(Y) = Yb com YA<c e Y qualquer, 
transforma-se em max g(Y) — Yb com Y(Al-— |) < (cl 0) 
donde YA<c e — YISO equivalentes a Y> O, portanto 
finalmente em: max q(Y) = Yb com YA<ce Y = O que 


é precisamente a forma dual do problema simétrico e 
porianto os problemas assimétricos e simétricos equi- 
valem-se. 


E com esta transformação que acabámos de ver, 
basta demonstrar o Teorema da Dualidade para os pro- 
blemas primal-dual assimétricos mas antes disso veja- 
mos certas relações entre o quadro inicial e o quadro 
final da resolução dum problema pelo Método do Sim- 
plex, nomeadamente no respeitante às matrizes corres- 
pondentes, às bases e às matrizes inversas: 
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a dn TD a TRA lis maiden tia CA e 2 Cit iii ati Cai ————— a 


INICIAL a 1 =3 O 
P 


PR. TIN Ss 0 
vector a€- P, 0W|)0-2 454 1 
Po QIO/O-4 30 


sair, P, 


0/0 =-1 3 0 
! 


* O 
P. P, — vectores coluna P; 


2 
O O — matriz A, base inicial (P,P,.P,) 
8 


—2 O — números Zk — Ck 


solução inicial vector a entrar, P, 


FINAL 

Base co Xo|X, X, X, X, Xj X, — vectores coluna X 

P, 1 4 Ef; 1 jô LF so é dr 0 

P, —-3 5l/, O 1%, *%, O — matriz X, base inicial (P,,P,P,) 

Ra O 1111 [8] O —3/, 10 1 ou B= 2 =1 0/0 
EO 0 5,250 a. 4 0 

solução final — 4 3 1 
os vectores iniciais P,, P., ..., P, iniciais estão expres- yYob= És Bo! b= Cc, X, = min fix) 


sos em termos da base final B=P,, P,, P, ou seja 
E = B KR,» como por exemplo se pode verificar para P., 
Pg = (2,0,8) = (Pa Par Pol (fe 2! 10) = (2,0,8) 

e igualmente se veria P, oub=BX,,istoé, A=B x 
ou B-'A=X e também B7!b=X,. A solução miíni- 
ma neste caso é c, x, =(1,-3,0) (4,5,11) = —11. 

E assim com a recapitulação destes pormenores 
e o esclarecimento de certas notações, podemos enfren- 
tar o teorema seguinte, 


3.2 — Teorema da Dualidade 


Enunciado: Se o problema prima! (ou dual) possui solu- 
ção óptima finita, então o problema dual [ou 
primal) também a tem e é min fix) = 
= max g(y). Mas se um dos problemas tem 
solução ilimitada, o outro não tem quaisquer 
soluções admissíveis. 


Demonstração: 


| a) Admitamos que o p.imal tem solução óptima 
mínima e obtida pelo Método do Simplex, isto é, 


min fix) = C, Xp & sendo mínima, por um dos teore- 


mas anteriores z — (cy; X — CC Xi — Cu 0 X, 
-Cl=(Z—C oleo l=0GA-—-ESO. 
Definamos agora Yo = (y0,y,0 ,...ym |) por 
Yo =€, B”!, Provaremos em primeiro lugar que é so- 
lução possível para o problema dual. Ora como X = 
=" A é Z=6,X-= 6 de Yo A< co Yo A- 
- c< O, vem: 
YoA-c=eBlA-c=e X-c<o 
ou Yo A-c< O E o valor da função objectivo dual, 
usando a definição de Yº e a igualdade Ka =B8"!b, é 
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Falta ver agora se Yº , além de solução admissível é 
também solução óptima do dual: 

Para quaisquer X e Y, olhando as formas primal e 
dual, vem: 


olv)=Yb=YAX<cX=t(X), de YASc 


ou seja g(Y) < f(X), desigualdade que vale também no 
caso particular: max g(Y) < min f(X), isto é, o máximo 
da solução dual é quando muito igual à solução primal. 
Ora já se viu que para a solução Yo do dual se tem a 
igualdade g(lYo)=Yob=c, X,=HX, ) = min f(X) 
logo temos 


max g(Y) = min f(X) 


b) Admitamos agora que o dual tem solução ópti- 
ma finita. Demonstra:emos que o primal também a terá. 
Para isso transformamos o problema dual para a forma 
primal e demonstraremos que o dual deste, com solução 
óptima finita pela alinea a) já demonstrada, é precisa- 
mente o primal original. Tem-se para o dual: 


max g(Y)=Yb com YA <ceY qualque: 


então max (Yb)=-min(— Yb) ou, — max (Y b) = 
= min (Y b) sujeita à condição YA + Y,I = c em que 
Y. é um conjunto de variáveis de folga não negativas 
e | matriz identidade. Ora Y é qualquer mas decompon- 
do em variáveis não negativas Y = y, — y,, vem a forma 
primal: 


min ((— Y, + Y,)b+ Y,0) sujeita a 


Y—-Y)A+rY,I=c com Y, > O para i=1,23 
e agora já se pode aplicar a alinea a). Falta apenas ver 
que o dual deste problema é o primal original. 


Multiplicando por — 1 a igualdade anterior vem 
—Y,A+Y,A—-Y,l=—coe para o dual: 
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mex (- cX)=—min (cX) sujeito a B7lg =B"! (Pisa Pg sssgiê 
(—- ALIAlL—I X <t— bIbIO) ou seja 


— max (— cX) = min (cX) sujeito a O cs Ri 
———————— 7? ja 
2h 
—AXS—b AX=b -1B 
| =ã A , isto é, B B =, 
| AX< o o PARTA com E 0 end 
| IX < O... Xma 


e este sublinhado é o problema original. ia ' 
tonde |, designa a matriz identidade com excepção da 


Il) Se o primal é ilimitado, então como g(Y) < coluna |). Então B BIB = Bl,, ou B = B|,, donde 
< fix), vem que g(Y) < — v, pois toda a solução dual Masai E (B| p= au E 
satisfazendo YA <c obrigaria a glY)= fIX) o que é = Ei : 
impossivel por f(X) ser Ilimitada inferiormente (de mim oh “por fp n SERIE der ealouiar ve Mundo 
f(X)). E da mesma maneira se vê que quando o dual 
tem uma solução infinita. o primal não tem soluções. 


=| !B7!. Mas decignan- 


e: 
1 Dm= das 
3.3 — Algoritmo do Método do Simplex Dual *k 
K. 
NOTA: Convém estudar em seguida o Método do Sim- O 1. eosimna = 
plex Dual pois além de resolver o problema dual X Ike Ba =E 
tem aplicações num algoritmo importante da E = pu ra 
Programação Linea: Inteira. Antes disso, porém, 1 B-!, o que equi- 
vejamos um facto importante para esse estudo: O O... 
A matriz inversa de cada nova base pode obter-se | X|k 
directamente da base anterior aplicando as fór- a 
mulas de eliminação completa, e aparece no O O... —-—& 
mesmo quadro no lugar dos vectores unitários ÉLk 
do quadro anterior. 
Seja então B= (P,... PP; qu Pu Ja base an- vala à aplicação das fórmulas de eliminação. De facto 
tiga e seja B=(P,..,P,,...P,)a nova base onde | designando por b,, o elemento genérico de B ! e 
P, substituiu P,. Ora por b;;o de B. 1, fazendo o produto E, BO 1, vem: 
las EA = 
8 B=B (PPP) = 
D,;j | D, 
1 EM parasse DO E, b. A para i +| :b, —  —— cem 
R Ly] x ,k E — ris 
O: Ô Ek IR 
“Jo o. ..o| evstoquex=B “P, Vejamos o exemplo seguinte que não é mais do 
que um extracto da resolução dum problema pelo Méto- 
O casi 1 do Simplex: 


B=[(P,P, P,), base inicia! 


elo ss q B*=iPytada 


- s 1 40 1 O 1/, 3 4 B=WafP) 
PL 3 0 1, o q ai 
Pp 1 0-4, 0-3, à à F'=mXGe=lo 3, 6 
Bo 6% mito cê O 


as igualdades são fáceis de verificar quer efectuando 
os produtos respectivos que; por aplicação das fórmulas 
de transformação, como aliás já se sabia pelo Método 
Simplex. Vê-se que a matriz inversa para cada base se 
vai fo:mando no lugar dos vectores unitários. 


3.3 — Algoritmo do Método do Simplex Dual 


Seja o problema primal: min cX com AX =b e 
X>0 


o seu problema dual será: max Yb com YA<ce 
Y qualque: 


Escolha-se a base inicial B=(P, P,,...,P, ) de vec- 
tores unitários de modo a que um elemento Xi o de Bb 


seja negativo, e portanto base não admissivel para o 


primal, e com c, X; ou Cc, E E, < c; para todo 
o j. Como se viu nô tedrema “da Dualidade ums solução 
admissível para o dual é: 
P=o Br! com gtry=c,B !b 

O objectivo é desenvolver um algoritmo para o 
problema dual de modo a obter a solução maximizante 
a portanto também a minimizante para o primal, como 
diz o teorema anterior. Para isso temos de ir escolhendo 
novas bases que dêem soluções cada vez maiores para 
a função objectivo até chegar ao seu máximo, se ele 
existi, e que além disso continuem a satisfazer as de- 
sigualdades YA<C, 


Designemos por B, as linhas de B- HA base não 
é admissível para o primal mas se o fosse diriamos que 
a sua solução era x,, = Bb (ou X, = B”!b) para 
[= TA 

Consideremos agora x, = min B, b< O que existe 

| 
por hipótese, e para os vectores que não estão na base, 
=p ;<c; ou 


4 =B, P, 8 su- 


isto é, aqueles que são tais que Cc, B 
For, <e, 


he Aa que pelo menos uma alisa X| 1 O 
Calculemos 4 da seguinte forma: 


calculamos todos os X| 


E AP 3 Ec it 
9 = min + ps Se O 


x X 
Xu<o bJ Eh 


Entrando o vector P,, seleccionado, para o lugar de 
P,, obtém-se nova base B que dará nova solução possíi- 
vel para o dual. B” 'bode obter-se como se viu pelas 


fórmulas de eliminação ou pelo produto EB”! Ba 
nova solução será 


E a a] ci a A 
Y=68,B = €, EB = 
= Es Ej= (Cj. Cy, Cj Cm! E, = 
B44 Ko k 
C,. C ] (— : - Ea l. Fei -+ 
|k *,k 
| Xm,k 
! Ea ko o m | Emo 
|,k |,k 
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acrescentando e tirando ao interior do parêntesis a par- 
cela necessária para termos 4; 


x 
Lk 
ac Cc, tc; vem: 
*Lk 
= E Cas IC; És 8), E 
= C,, Cu, oem Ep een E — 8 (0,0 1 ,+0) = 


= Fac 9 B, por ser de vectores unitários a base inicial. 
ou seja 


Y=Y,-08B, 


] 


sendo o valor da função objectivo, por ser B, b=«x, o 


Claro que a nova solução do primal é X=B !b 
e que também será óptima se todas as componentes 
forem não negativas, X> O. Se isso não suceder é por- 


que pelo menos uma componente Kyo E B, b<o e 
nesse caso repetimos o trabalho anterior determinando 
o vector que entra e o que sai para gerar a nova base 
atô à solução do dual e consequentemente do primal. 
Estas soluções podem não existir, se for ilimitada a 
solução do dual, e isto ocorre se forem todas não nega- 


tivas as 5 = B, P, podendo agora escolher-se um 


à > O qualquer. De facto, sendo X o =. Era e, 
cresceria tanto quanto se quisesse. 


1,0 


Capítulo 4 
PARTE PRÁTICA 
4.1 — Esquematização dos Algoritmos 


a) Método do Simplex. Recordando aqui o quadro 
inicial exposto atrás, podemos esquematizar este algo- 
ritmo da seguinte maneira; 


1.º Apresentado o problema, transformamos o 
sistema das restrições num sistema de igualdades 
usando variáveis de folga e de modo a serem não nega- 
tivos os termos independentes, P, , no Quadro Inicial. 
Se o objectivo for maximizar, mudamos para minimizar 
com max f=—- min (— 1). 


2º Se conhecermos um conjunto de vectores 
linearmente independentes tomamo-los para base inicial, 
se não, tomamos para base o conjunto dos m vectores 
unitários (P,, P., .... P,, ) formando a matriz identidade 
tm X m), base admissivel, que se pode sempre arranjar 
com a utilização de variáveis artificiais e construímos 
o quadro inicial. 


3.º Assim temos a primeira solução inicial X, = 


ER Ra | 
HE Po ÀXio X50º «+ Xmo) com x, >0, vector 


ti 

: =| = 
3 Cc; Xijo * ainda X; =B P, = 
i—1 


coluna, 
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= (je Xy jr == Xm |) vectores colunas também. 


4. Calculamos os elementos da m + 1 — ésima 


mm 
linha do quadro, z,— ce = > C; Xjj —€, para 
i=1 | 
j=1,2.....n e colocamo-los nas respectivas colunas 
| Os z, — c, referentes aos vectores que estão na 


bases são sempre nulos. 


5.º Se os Z; —€; são todos não positivos então 
X, 8 Z, é a solução óptima e paramos. Se não forem 


passamos ao número seguinte, 


6.º Consideiamos os númecos Z, —€; >0 e 
para diminuir o número de situações necessárias, con- 
sideramos max, (z; C; ) = Ze — CG Se aparece- 
rem indecisões tomamos o vector coluna de menor in- 
dice j. O vector P, é o vector a entrar na base. 


7º Do vector P, analisamos as suas componen- 
tas Xik Se forem todas Xi k <O o valor da função ob- 
jectivo torna-se ilimitado e paramos. Se não, das 
X; 20 calculamos 


hj Xo 


ah 


5 A 
e o vector P, sai da base para dar lugar a P,. 


8. Passamos a construir o quadro seguinte onde 
a base apresenta o vector É” no luga: de P, . O ele- 
mento pivô é o X| y Cruzamento da coluna k com a i- 
nha |. Todas as linhas incluindo a m + 1 — ésima são 
transformadas do primeiro quadro pelas fórmulas de 
eliminação, o que equivale a exprimir todos os vectores 


P, que não periencem à base em termos da nova base. 


9,º Depois de elaborado o segundo quadro temos 
a solução admissível X', = eq: js X ko 205 ma) 
de 


Com X;o De 2,=%”- -X4 )8 passamos 


Xik 


ao número 5 para continuar o processo, 


NOTA: Se o p.oblema contém variáveis artificais, em 
número de m ou menos, sabemos que a função 
objectivo será min c x, +e X + ... + 
E CÃ£o TWXds FW go Tvs EWX 4a 
e a base, artificial portanto, será (P,.ge 
é | para o sistema aumentado. Temos 


nr mm 
Z.—C=W Wo «o W).P; —G, = 


J 
na 
= W > Xj;— C com uma parte em W e ou- 
i=1 
tra sem W. 


Assim fazemos uma alteração ao quadro já descrito 
criando a m + 2 — ésima linha, colocando nela os coe- 
ficientes das partes em W e na linha m + 1 as partes 
sem W, sofrendo também aquela linha as t:ansforma- 
ções habituais de quadro para quadro, 
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A selecção do vector a entrar na base faz-se a 


m 
partir da m + 2 — ésima linha e toma-se o max 5) Xi 
j=1 

como vulgarmente e segue-so o algoritmo descrito até 
que todos os vectores artificiais saiam da base ou até 
que não haja nenhum elemento positivo na linha m + 2. 
No primei.o caso a linha m + 2 tem todos os elementos 
nulos e a base obtida é admissível para o problema ori- 
ginal e portanto seguimos agora o Método do Simplex 
vulgar. No segundo caso há duas possibilidades, ou o 
elemento (m + 2,0), isto é, a parte artificial do valor 
da função objectivo, é maio* que zero ou igual a zero. 
Se suceder a segunda, (m + 2,0) = O, estamos em pre- 
sença duma solução admissível com pelo menos um 
vector artificial, degenerada portanto, e os cálculos con- 
tinuam. Seleccionamos o maior dos elemetnos da linha 
m + 1 que ficam na mesma coluna dos zeros da linha 
m + 2 para escolher o vector a entrar na base. Atinge-se 
a solução minima quando não houver mais nenhum ele- 
mento positivo por cima de zeros da linha m + 2. Em 
ambos os casos os elementos da linha m + 2 não são 
t.ansformados de quadro paraquadro. Se a solução final 
contiver variáveis artificiais elas serão nulas. 

Se apenas tivermos necessidade de algumas variá- 
veis artificiais e não das m referidas, o processo cor.e 
da mesma maneira, 


b) Método do Simplex Dual, 


Valendo-nos do conhecimento do Método do Sim- 
plex, torna-se fácil esquematizar o procedimento do 
Método do Simplex Dual pois eles assemelham-se mui- 
to nas suas linhas gerais. 


1, Escrevemos o problema na forma: Min cX su- 
jeito a AX = b com X5>0, com o seu dual Max Yb 
sujeito a YA<c. 


2.º Tomamos uma base inicial (P,, Pa... Ph) 
tal que pelo menos um elemento B”!b seja negativo e 
com €c, — and Ps 4 c, para todo o j, formando o qua- 
dro inicial. Temos a solução Ya = és B— 1 eo valor 


da função objectivo c BB 'b. 


3." Sendo B, as linhas de B- pá calculamos 
= min B b< O e é eliminado da base o vector P, 
i 


Xo 


4." Determianmos as componentes X|j < O per- 
tencentes a vectores colunas com m + 1 — ésimo ele- 
mento tal que YO PR < c; e passamos ao número se- 
guinte. Se forem todas x; =» O o processo acaba e não 
há solução finita. 


= min 
ALi<o Li 
8 o vector P entra paia o lugar de P, 


&.o Calculamos 


6.º Dispondo de nova base construímos um novo 
quadro em que os vectores não pertencentes à base es- 
tão expressos em termos da nova base, usando para 
tal as fórmulas de eliminação completa, 
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7º Se ainda existir elguma componente Kia < 0, 
passamos ao número 3.º e repetimos o ciclo, Se for 
Xj02 O para todo o i atingimos a solução óptima para 
o primal constituída por essas mesmas componentes e 
portanto também solução dual que é C, B'! 


4.2 — Problemas Resolvidos 


1º Minimize x, — 3x, + 2x, sujeita a x, > 0 e 


TSM * 2% =7 


— de e AX PR = 12 


— 4x, + 3x, + 8x +x,=10 


imediatamente se vê que dispomos dos vecto es unitá- 


rios P, P, P, em número de m (3) para formarmos 
a base inicial e podemos formar o primeiro quadro, 


=, 0 2 DB 


PRF 


3 P, O 10 O —4 3 0 8 1 
4 O 0 = RR a O 


em que c, =0,=€C=0,X, = (M XX) = (7,12.10) 
2, =>cX, =0 e a linha 4 é constituída pelos ele- 


Ea 


mentos Z, — c, = cr P, -C, para | E Lcd sea 8. 


Max lz, — c, )=2z,-— c = 3>0 eportanto entra pa- 
J 

“E : O y 12 

raa baseo vector P,.6,) min) — = min ETR 

Xxis>0 Xo 4 


10 12 
| 


má e sai portanto P,. 


Passamos ao quadro seguinte aplicando as fórmu- 
las de eliminação 


2 PF; TS 4d P-11, 1 44, 0 O 
3 Po O 1 O =8/, 0 —3/, 1 
4 2 O %a GS 7 O 
e a nova solução é Ka = (Xj Xi Xe) = (10,83, 1.2", = 


| 10 
= — 9. Max (z;—6,) o 50. 6 ma! Po =" 


e sai P,. Transformando novamente este quadro, vém 
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O O —=4/,-12/, 0 


A nova solução é X' , =(x,, xy x) =(4,5,11) e 
Eca = — 11. Como max tz, = 6; |) = O esta solução 
é a minima. 


NOTA: Vê-se que são todos nulos os Z;— €, para 


ji=2,3, 6 mas se algum z;—6,=0 para 


j-H2, 3, 6, então um desses vectores podia 
entrar para a base e obtinha-se uma nova solu- 
ção minima mas sem alterar o valor da função 
objectivo. A solução final seria a combinação 
convexa dessas soluções, 


2.º Maximize x, + 2x, +3x,— x, sujeita a 
X; > 0 e 
à FEM TSM, = 15 
ds dt 2 + Dx; = 20 
Mr 4 TS TO 


em primeiro lugar transformamos o problema para mi- 
nimização: — min —- x, — 2x, — 3x, + x, Como dis- 
pomos apenas do vector unitário P, necessitamos de 
dois vectores artificiais, isto é: 


— Min —x,— 2x, —- 3x, +x+tw+w com 


X 2% IX + w = 15 
ER TE a TEM + w = 20 
Her Mgmt; = 10 


e formamos o primeiro quadro: 


Co.-l-L=-3 7Tlw w 


| Base c 
TOP, w 
A Pa W 
3 P, 1 
4 10 2 4 4 010 0 
5 sb 33 80140 0 


Notar os elementos da linha 6 (m + 2), coeficientes de 


5 
w, como por exemplo z;, — c, = w b2 XxX —6G>=>wWw+ 
d=10 0 


25 


+liw+r1i-—-(l—1])=2+3w:2 na linha 4 6 3 na 
linha 5. A solução inicial é X, = (x, Xp x) =(15, 
20, 10) e Z, =10+35w. Max im + 2,j) = (m + 
+ 2,3) — 8 e entra o vector P,, 6 = 20/5 (valor mi- 
nimo) e portanto sai P,. Aplicando as fórmulas de eli- 
minação, vém:; 


| Base c Po FP Fr fa Fi wu 
1 PP, w Bt O 6 1 
> Po =3 à 4% 4 6 0 
2 IR. 1 Ds EC O 1 O 
4 = fd UU ! O 
5 3 2h tt D O O 
em que X, = (x; X% x)=13,4,6)02Z2,=—-6+ 


+ 3w. O máximo da linha m + 2 é 7/5, entra P,; 
9, =3/"/p o sai P,: 


Pp,  —2 bj, if 1 0 0 
Pp, —3 v/, 4, 0 1 0 
P, 1 sf,  % O O 4 
—s0 4 0 O 0 

0 o o 0 


Como (m+2,))<O0O para todoo je (m+2,0)=0 
temos uma solução possivel para o problema original 
que é X = (x, xy x,) = (15/7, 25/7, 15/7) e Z/ = 
= — 90/7. Como já não há nenhum vector artificial na 
base, passamos a aplicar apenas o Método do Simplex 
para conseguir a solução óptima. Entra P, e sai P,: 


FP see ta O 9 (6) La 
ur =:3 o Oo O 1 a 
F, a to 1 O O Ela 
=1% 0 0 0) = 
E como z, G < O para todo o j e apenas t; — 
= = O para |= 2,3, 1, obtemos já a solução ópti- 
ma que também é única. X" =x, xy x) = (5/2, 
5/2, 6/2) e Z",=— (— 15) = 15. 
3.º Minimize  — xp + 150 x, — 
1 
—k Tt Ox, sujsttoa x, > O e 
50 
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fa K4 =60% "— ask Tt RT Es ms 
at — 80% Yet, FIA, + x = 


Ent 


+ oo 


Ka + x = 


a base inicial é evidentemente (P., P,, P.) e formamos 
o quadro: 


C=»=4/, 150-1/,, 6 O O O 


Rê] | a 


| Base ce P P Es Po Pe Pa mr P; 


1 PL 0/0 1»,-60-14,, 9 1 0 0 
à Pp OI O 1/, —90 — ja 8 0 4 q 
2 PO DA1moO Ga 
4 O 3/,-150 1/4, —68 O O O 


Logo se vê que entra P, (z, —c, = 3/4 é o máximo) e 
su.ge egora indecisão quanto ao vector a sair, P, ou P,? 


1 1 
Usando a técnica descrita em 2,3., min ( 
4 
1 1 
] 2 2 nada resolve mas indo a |j= 2, min 
1 1 
| — 60: :—90: —) =Xya/X = — 240, logo 
4 2 | 
sai O primeiro vector ou seja P 
, PF, ="Ê; O 17—-240 —3/ 36 4 0 O 
* E É) O O o “ho cta-2o 10 
3 P, O 1 O O 1 O 001 
4 O O SO “fas -=3 O O 


aqui entra P, (30 é o máximo) e sai P,: 


1 P,-47,| 01 0 54, —-84 -12 8 O 
2 P, 150 | 0 017/40 
3 P. O 100 1 o 0 01 
4 O 0 0 “4,18 -1-1 0 


temos novamente dúvida pois x,., = X,y = O mas des- 
fazemo-la calculando min (1/5/,,:0/1/.)]) =0 o que 
faz com que saia P, 


TO PL —a3p O 1-160 0 —4 44, 8/4 O 
2 P, -!/u |O O 500 1 —250-100/,%97, O 
3 P, 0 1 O —500 O 2501007, —50/, O 
q 0 0 -40 0 2 Syj —It 6 
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também aqui não há qualquer dúvida, entra P, e sai P,. 


2 


| 2 
1 P, —3), ET 1-168 O O —ty, ci T 
à Elio UR: 0 1 o O o 1 

1 1 
Ss P, 5 | SEO o => q q Rea 
4 Rai d O -96 0 6 tj -uy Dra 

125 : 125 


entra P, e sai P,: 


2 PF, Ficá 1 O O 1 o o õ0O 1 

| 3 3 
3 P, (6) 100 “=15 BWj, V=f; TOS 
4 “20 O 15 O 0 na “20 
Todos os elementos z —€6 são negativos e nenhum 


dos respeitantes a vectores não pertencentes à base é 
nulo, portanto a solução encontrada é óptima e única: 
X = (Xe Xp X) = (1/15, 1, 3/100) a 2z= — 1/20. 


4.º Resolver usando o Método do Simplex Dual: 
Max 7x, — 7x, + 2x, +x, + 6x, sujeitoa x, > 0 e 


So mta ts, =—3 
Zum Xi + rm; =" 4 
=x FD - + x = 12 


Para termos o problema na forma canónica temos de 
transformar em minimização: — min — 7x, + 7x, — 
—2x,—x,— 6x, e a base que se deve escolher é 
B=(P, P, P,)=B1 pois existe um elemento de 
P, — B'! b negativo e todos os C, B-1 P<c,2Z; <€; 
ouUz; — €; « O, como se vê no quadro seguinte. Claro 
que a base não é admissível para o primal: 


G ip EE =] =B iu 


Base c Fa PB Pa PP Pk ER, 
4 Pr SJ 4 = 414 O 0 
2 Fr =B 4 2 1 O a 1 O) 
s Pi O pa om Ss 0-3 0 1 
4 [8 =11 0 =08 0º =1 o O 


A solução do dual éY, = c, B1i=(—-2, —6,0).B1= 
=(-2, -6, 0).Min B, b=x,,= —3, sai o vector 
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P,. Paraj £3,5,6,x; =B,P; =(1,0,0). (P, P,, 
P,) ei (3, =], — 2) * (—, = | — ), 
9 = min dos mmij= —— | == 
Xij>o a Ma = 
1 1 alo 
ns O, = ———— logo entra P,: 
a 2 dci , 


7, ic D/, = da O OQ 1 


Como já não há nenhuma x, , positiva e z; — c, são 
todos negativos, chegou-se à solução óptima para o 
primal: x, = 3/2, x, = 5/2, x, = 33/2 e com o valor 


da função objectivo — (—33/2) = 33/2. 


A solução óptima para o dual será Y', =c', B!= 


-1/, 0 0 
=(—1,-6,0) 1/1, 1 0 = (—5/2, +6,0) 
-81. 0 1 


ou seja: 
y=-5/2,.y=-6,y,-0 


O valor da função objectivo do dual é g(Y) — 
e Bib=Yb=(-5/,=6,0)./ -3)=- 


[a 
2 


5.º Resolução do problema do Ourives 
Maximize Z=3x, + 4x,  sujeitoa x, x, >0 é 
x, + 2x, < 28 
3x, + 2x, « 52 


A forma canónica será: — min(— 3x, — 4x.) sujeito 
à Ko + 08 


a assim: 
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Po —4 Ss qq 4 MH, 
P. —3 12 1 O =27, djs 
-98 0 O -3/, df, 
Logo x,=8, x,=12 e Max z= -— (— 98) = 98. 
PARTE II 
Capítulo 5 
PROGRAMAÇÃO LINEAR INTEIRA 


5.1 — Introdução 


A Programação Linear Inteira (P. L. |.) estuda pro- 
blemas semelhantes aos de Programação Linear Geral 
mas em que se exige que todas as variáveis envolvidas 
tomem exclusivamente valores inteiros. Muito embora 
haja um outro tipo de problemas de Programação Linear 
Inteira, o chamado problema misto, em que uma parte 
das variáveis têm de ser obrigatoriamente inteiras en- 
quanto as restantes podem tomar valores contínuos, será 
nosso objectivo dedicar-nos aqui ao estudo do primei:o 
tipo de problemas. 

E é grande a aplicação prática da P. L. |. pois que 
se há problemas onde a Programação Linear geral nos 
dá como valor da ce.ta variável o número 120,3 (por 
exemplo hectares a semear ou número de trabalhadores 
a utiliza.) onde 0,3 pode ser desprezado, são também 
vulgares as situações em que se procura por exemplo 
o número de navios a construir ou a desempenharem 
uma certa tarefa e nos aparece um valor não inteiro 
cuja parte decimal é importante demais para ser des- 
prezada, 

Claro que com determinados artifícios nos dados, 
estes problemas poder-se-iam resolver pelo Método do 
Simplex num número finito de passos mas infelizmen-e 
o tempo de resolução tornar-se-ia demasiado longo ou 
os resultados obtidos seriam de pouca confiança em 
virtuda dos arredondamentos necessários. Ássim desen- 
volveram-se e continuam-se a desenvolver muitos algo- 
rimos que oferecem possibilidades em variadas aplica- 
ções sempre c.escentes, conquanto nenhum deles tenha 
a gensralidade do Método do Simplex. 


Uma vantagem de que goza a P. L. |. é poder-se 
introduzir no modelo do problema relações lógicas com 
a utilização de variáveis booleanas, 


Neste trabalho serão apresentados, com a respec- 
tiva programação em linguagem A. P. L., apenas dois 
algoritmos que nos pareceram de mais generalidade e 
alcance prático: o método de cortes, algoritmo Gomory, 
para variáveis inteiras quaisquer e um método de Enu- 
meração Implicita, algoritmo Balas onde as variáveis in- 
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teiras são obrigadas a tomar exclusivamente os valores 
O ou 1 e que serão sempre representadas por 5:0< 
< 5 < 1 com jôiinteiro, variáveis booleanas portanto. 

Além disso com a transformação de quaisquer va- 
riáveis inteiras em vaiiáveis booleanas ampliamos a ca- 
pacidade do algoritmo de Balas, pois nesse caso passa 
a poder resolver todos os problemas resolúveis pelo 
algoriimo da Gomory, tendo apenas como limitação o 
número exageradamente grande de variáveis boolenas 
qua por vezes aparecem naquela t;ansformação. 


5.2 — Aplicações e Modelos 


Logo no início do desenvolvimento da Programação 
Linea: muitos problemas foram enunciados em termos 
de Programação Linear Inteira por Dantzig, Hirsch, Tu- 
cker, Wagner e out:os. No entanto foram os trabalhos 
recentes de Gomory e Balinsky, que apresentou uma 
grande variedade de métodos, Benders e muitos outros 
que abriram novas vias à resolução prática de muitos 
problemas de P. L. |. 

Entre os p:oblemas mais frequentes que não po- 
dem ser resolvidos satisfatoriamente pela Programação 
Linear geral em virtude dos inconvenientes já aponta- 
dos, contam-se os problemas de organização de trajec- 
tos terrestres, maritimos ou aéreos, de sincronização 
do semáforos, de carga fixa, de programação de ma- 
quinaria, de stoques, do caixei:o viajante, etc. Mas mui- 
tos outros existem totalmente fora do alcance do Mé- 
todo do Simplex, nomeadamente problemas em que o 
conjunto das restrições é formado po; sub-conjuntos 
de restrições mutuamente exclusivas (condição «our), 
já que teriamos de efectuar a união dos conjuntos con- 
vexos das soluções de cada sub-conjunto que em geral 
não é convexa. No entanto é precisamente com a in- 
trodução das condições do tipo «ou» que se tornam 
resolúveis certos tipos de p:oblemas onde a função a 
optimizar ou algumas das suas restrições não são li- 
neares, 

Não é possível em P. L. |. formar um modelo ma- 
temático único englobando os diferentes tipos de pro- 
blemas, contudo para os resolúveis pelo algoriimo de 
Gomory esse modelo só difere do conhecido para a 
P.ogramação Linea; ge:al pela restrição adicional que 
obriga as variáveis a tomarem valores inteiros: 


min HUX)=c x +..+e X, 


com 


dm1 


X Fo + à mn An Dm 


e todas as componentes de X inteiras e não negativas. 
ou mais sintaticamente 


min cX sujeito a AX =b e X > O de componentes 
inteiras. 

Para os problemas de variáveis booleanas considera-se 
o seguinte modelo: 
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min cX sujeito a AX <b, 0O<x, <17, x inteiro e 


com c = O 

NOTA: Para uma melho; compreensão de certos aspec- 
tos dos assuntos expostos, especialmente no 
respeitante a notações convém recordar da PAR- 
TE | deste trabalho tudo aquilo que possa 
contribuir para qualquer esclarecmiento, já que 
as matérias se encontram interligadas. Assim 
só será feita referência às ideias e notações 
que surjam aqui pela primeira vez. 


Capítulo 6 


ALGORITMO DE GOMORY E SUA PROGRAMAÇÃO 
EM LINGUAGEM A.P.L. 


6.1 —— Generalidades 


No método do Simplex começa-se por encontrar 
uma solução de ponto extremo e ao ir seleccionando 
novas bases encontram-se novas soluções de ponto ex- 
tremo até encontrar aquela que optimize a função ob- 
jectivo. Ora o nosso problema consiste em encontrar a 
solução admissível de componentes inteiras que mais 
se aproxime da solução, geralmente não inteira dada 
pelo Método do Simplex. 


Uma primeira tentativa de resolver o problema e 
devida a Gomory foi acrescentar restrições ao modelo 
do p:oblema represntando cortes no conjunto convexo 
das soluções, Acrescentando estas restrições ao modelo, 
o problema conserva a sua possibilidade de optimiza- 
ção mas nem sempre a sua resolubilidade pelo Método 
do Simplex. Mas aplicando o algoritmo do Método do 
Simplex Dual a resolubilidade aparece e se a solução 
óptima ainda não é inteira introduz-se uma nova res- 
trição e aplica-se outra vez o algoritmo dual até se che- 
gar à solução inteira ou à certeza dela não existir. 

Mas este processo foi ultrapassado por um outro 
mais eficiente e também devido a Gomory em que se 
parte dum problema com dados inteiros, sempre possível 
de consegui:, e que por aplicação do Método do Sim- 
plex Dual sujeito a certas modificações, todos os qua- 
dros conservam a natureza inteira dos dados até à solu- 
ção final. Em linhas gerais este processo consiste no 
seguinte: 

Começa-se por construir o quadro inicial do Sim- 
plex com todos os dados inteiros e escolhe-se uma 
solução admissível para o dual mas não admissível para 
o primal. Para isso toma-se uma base B = (P,, P.. ..., 
P) nas condições exigidas pelo Método do Simplex 
Dual, isto é, tal que pelo menos um elemento x, , = 
= Bb seja negativo e todos Z; —€; < O, sendo 
então Y,=co Bb, respectivamente, a solução e o 
valor da função objectivo do dual. Em seguida acres- 
centam-se às restrições do modelo mais as seguintes 
aquações: 


x, =x; para j=m+H1, 


«o N comx, =0ec,=c 


J 
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Esolhemos agora uma solução possivel que contenha 
n variáveis não negativos em vez das m habituais para 
o Método do Simplex e que por não serem linearmente 
independentes nos dão uma solução interior ao con- 
junto convexo das soluções do modelo inicial. 

O algoritmo desenvolve-se agora procurando um 
processo que possibilite passar deste ponto interior para 
o ponto óptimo de coordenadas inteiras. Assim intro- 
duzem-se novas variáveis à solução de modo a conser- 
var a natureza inteira dos quadros. Isso consegue-se 
definindo certos planos, restrições em termos das va- 
riáveis não-básicas e dum vector coluna de folga Pas! 
vector este que além de manter inteiros os dados ao 
longo dos quadros, vai permitir deslocar o ponto inte- 
rior inicial para mais perto da solução óptima inteira. 
E com a aplicação sucessiva destas restrições se chega 
à solução óptima inteira, 


6.2 — Escolha do plano de corte 


NOTA: Cada elemento x,, dos quadros, inteiro ou não, 
pode ser escrito como múltiplo dum inteiro 
mais uma parte fraccionária chamada resíduo: 


% - 
= E 4 |. a y 


q 
=) == xd! 
A 


em que O < prt E Ns Dj; é inteiro, r,,;6 o resíduo, 
à é um número positivo que se irá determinar e [...] 
significa «o maior número inteiro contido em». Evidente- 
menta que se X,,;/A< a sua parte inteira é b, ,< 0 e se 
AD 1, [1/A] = 0. Para melhor compreensão vejamos 
os exemplos: 
dps E 3,2 e seleccionou-se A= 0,5, então 3,2 = 
[3,2/0,5] «0,5 + "1; OU 3,2 = [6,4]* 0,5 + 4) 
donde 3,2=6xX0,5+r,, > n,;=0,2 e portanto 
como era de esperar: 3,2 = 6X 0,5 + 0,2 e 0,2 < 0,5. 
XE — 26 e h=0,2, então -Z5= [- 2,5/0,2] 
0,2 +14;= [—12,5] +0,2 + 1, ;= —13 X 0,2 + 1; 
= = "26 + 1) donde 17 = 0,1 < à como o exigido, 
Tomemos então uma equação com alguma compo- 
nente X, <O Seja a | — ésima: 


: Et E 


dio TO Pas Tlm o pmti  Hm+1 

-+ : 

E qe FM Ri (2.1) 
com Xm1 Xm+o! ---« XpVariáveis não-básicas e x, básica. 


básica, Aplicando a escrita exposta na nota a equação 
anterior vém: 
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ou ainda: 


IE 


igualdade esta que se verifica para qualquer valor inteiro 
positivo das variáveis como na | — ésima equação inicial 
e onde o segundo membro é positivo por se:em positivos 
os resíduos. Designamos o conteúdo do parêniesis () 
po: x, 4 « (variável correspondente à nova coluna P 
a com s a designar o número de planos de corte): 


E Xo | ha *j,j 
sara | : | = E ir ms pas X; (2.2) 


k um 
que se vê ser inteira pelas condições postas e também 
não negativa pois se o fosse, como r,, e À são positivos 
viria o segundo membro da penúltima igualdade negativo. 
Ora é precisamente este inteiro não negativo x, |. 
que se introduz como nova variável. Se restringirmos 
à > 1, como [1/A] = 0, vém para x, ..: 


n+ts 


| É base 


XL *m+1 | *n 
x e aq X — q — x 
n+s , N m+ 1 N ! E 
au, pela definição de b, 
Xy+g = Dog Db, m+1 Amy Te D, X, 
e finalmente 
Bjo ED mes + dC É Dia a E uag com A > 1 
(2,3) 


que é a restrição usada como plano de corte e que é 
satisfeita para qualquer solução inteira visto que qual- 
quer bj; e x,4+s São inteiros. Além disso como se 
supôs é X o <<, bjo < O e também alguma variável 
não básica X|j é negativa para que o Método do Sim- 
plex Dual possa ser aplicado pois se não o problema 
não seria resolúvel, 

Se à fôr suficientemente grande, qualquer Lx, ; / 
fh] para X,j < O dará b,; = — 1 e então a equação 
de corte, agregada na linha n+1 do quadro inicial dis- 
porá dum elemento — 1 que pode ser usado como pivô 
Parece então fácil escolher À mas sucede que quanto 
mais pequeno fôr À maior será o aumento do valor da 
função objectivo e portanto maior a rapidez com que 
se obtém a solução óptima (estamos a usar o Método 
do Simplex Dual e portanto 


£ E £ — 
p i k | 
& = min == Ea = >» O 
xico Xi Xi od 
E Tb=Y;b kg) 


De facto esse valor é agora 
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E jp = bo=x + 
"o. us: 1 Dj o o X 9.0 + X ok Lo * 5.0 
x 
RE 
+ iz, — 6] 
Iá [o h 
o que prova a afirmação feita e onde x, ,, é o valor 


inicial da função objectivo e P, , o vector coluna pivô, 
deve ser escolhido do mesmo modo que no Método 
do Simplex Dual e de tai modo que 


. £| = “(ga : En Xo.k E |. E Cj 
= min —— = min E a e ai 
aljoo X|j xlijco Ri 


6.35 — Selecção de » 


Como já se viu há conveniência no aparecimento 
do elemento pivô —1, cruzamento da linha restrição de 
corte com a coluna k, isto é —1 = Dj. = [x 1 / 4] 
de modo que pela definição de 5: 


A 2 o 7 
[xp /A] = Ex; /A] Dj 
õ it 
para x; < ou EE € b,; 


Designando ago.a por m, o maior inteiro para o qual 
se mantenha 


virá: 


e sendo x, < 0, o à mais pequeno a satisfazer esta 
desigualdade e que simultaneamente faça pivô a coluna 


P, ao conier o elemen'o b,, = — 1 será para cada |: 


Mas para que efectivamente apareça o elemento pivô 
= E Ú A; escolhido terá de ser pelo menos igual ao 
maio: dos o isto é: 


À = max A, 
alje o 


NOTA: Se o cálculo de à conduzir a A=1, isto é, 
À Ed = = Xiy = 1, temos de arranjar uma 
restrição de corte que não altere a base, ou me- 
lhor, que não a destrua. Para isso ela terá que 
ser função unicamente das variáveis não básicas 
o que só por si implica que À > 1 já que terá 
de se anular x, em (2.2). 

Sabemos que À não tem que ser inteiro e 
que para P, temos mp, AM =X pe 
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À > hp O que faz com que Ep = Lx, /h] = 
= — 1 como é necessário. 


6.4 — Esquematização do Algoritmo 


1.º Toma-se uma sclução admissivel dual e for- 
ma-se o quadro inicial com elementos inteiros e juntan- 


do as n — m restrições x; =x paraj=m+1,...n 


com x; >0e c'; = €; 
2º Toma-se uma equação (a | — ésima) que te- 
nha x,, <0Oetal que x ,=min x;,. Se não existir 


| 
nenhuma componente x, < O então a solução dual é 


já a óptima e inteira. 


| 


fm 
do que 


Calcula-se o vector coluna P, , pivô, de mo- 


— (Z,. -€, ) = mir 
aljeo 


[8-5 | para j>1 


Se não existirem X|; <Oop-oblemanão é possível 


7 


4.º Para cada j onde x,, < O, calcula-se o maior 
1 


inteiro m, tal que 
— lg —6 
— 4 — Sp 1 “To 
k k m; 
9.º Define-se À, = — (x,;/m, ) para todas as 


6.º Determina-se À = max A Se max A; = 4 


| 
toma-se um A > 171 qualquer mas tendo em menta as 


considerações sobre a grandeza de A vistas no núme- 
ro (2.3). 


7º Costroi-se a restrição de corte (2.3) e junta-se 
ao quadro inicial na n+ 1 — ésima linha, ficando a 
nova variável Xoys à corresponder à coluna P,.. 
(onda s representa o número de restrições que se vão 


construindo. 


8." Aplica-se a transformação do Simplex por in- 
termédio das fó-mulas de eliminação completa incluindo 
a linha n+ 1 e tomando Dj k = — 1 como elemento 
pivô. Assim fica introduzida a nova variável não básica 
Xn+s * como o elemento pivô era — 1, o quadro se- 
guinte será inteiro como o inicial e aplica-se novamente 
todo o processo, 


6.5 — Exemplo resolvido 


Minimize 10x, + 14x, + 21x, sujeita a x >0 e 


dx + tia rox, = 12 
mt 2x + 7x, Ra = 14 
8x + 6x +3x — x, = 10 


Todos os elementos são já inteiros; multiplicamos cada 
equação por (— 1) para obtermos termos independentes 
negativos como exigido no Método do Simplex Dual; 
tomamos para base P,, P,, P, admissível para o dual 
e inadmissível para o primal; acrescentamos as equa- 
ções X'j =x; para j=1,2,3, construímos o quadro 
com mais a linha n + 1 =7 e a coluna P, para serem 
p cenchidas com a restiição de corte e com a nova 
variável: 


c 10 14 21 O O o 10 14 21 
i Base c Po Po oP Pa Pe Pk RP PP Pal” 
a)t O O —TO —14 —-21 O O Õ O O O à] 
| og --d 1 6 06 d dO oleo 
=]. —% = ij O 1 ê ê O O O 
-10 -9 -6-3 0 60 1 6“ » olo 
O = O [6] O 6) O 1 8) O O 
b) O O = o O O O É) 1 o |0 
o o 0-1 o o o o o1“Io 
-4 1 1-2 0 0 oo “OO | 1 
onde a) é a linha dos Z —6G 8 b) conjunto das res- — 21 
rições: —x, +x,=0;-x +x,=00 —x, +x,= ET e Ps nho 
=0. | 
De P, toma-se a mínima componente, —14, e seleccio- é Cy z & 4 ou 16€ 10 14 l 214 
na-se a 2.º equação para gerar a restrição de corte. e ” qm m om MA 
Então dono 
——— = min "> 4 = min REAR ca im, Ma m,) < iq comme o: (DD e O 
= 1 pin ( -— ) X = 1 — 1 : É a a = 10 ' 10 , 10 | 
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= Xi Pa É Ts 7 ou —4&=-x —x,— 2x, +x, que se colocou na 
= max À, = max —— = max | 4 — 4 — )|= —— ' | : À 
m, (1 OBA É linha 7 do quadro anterio.: x, 4 = x,. Aplicando agora 
Logo a restrição de corte é: as fórmulas de eliminação X | Mi = (x; XD 
A us = 2 ni SA = ij xp para j=0,1,....7 e para cada i=0,71,.... 6 e 
e ed = pp == IM tl = Ka SE, a A 
7/2 7/2 T/2 1/2 . X pj EX Xp ac.escentando a coluna P,, vém: 
Base | Po P, P, P, E Pp, Ps P 1 Pa ZA P. P, 
O 40 D 4 1 O O O (6) O O-10 O 
t P. O 2 O —3 Z 1 O (é) O) o B]=8 8 
2 P. G|-ê O Do = í8 1 O O O 0) -2 0 
S Pr O 25 O 3 15 O O 1 O O O]-B 
4 P, 10] 4 [ô) 1 2 O O O 1 GD Bl] [ô) 
5 P', 14 D =1 o O O 6) O 1 O o oO 
E PP, 4 O O + 8) O O O O 1 O OQ 
F; 4 O O —1 O O 6) O O 014 1 
Temos x,)— —6B e portanto a equação geradora é a E = E | = | 
segunda, e determina-se a coluna pivô: és min ( Fai -] = 9 coluna pivô P, 
é PA O A | 5 24 
-4 = min ( -. — ) =1 5 Po to < (mo) Sm =1 e poranto = 7 = 24. Lo- 
| Er E 1 10 s — 4 3 
coluna pivô Im, m,) < == 4 = (1,10) 5 go restrição de corte será =) —— [ix 
24 24 j 
3 P. 
ANA Ls Ê , E S — 24 | 15 
À = max a ' ET = 3 Então a nova restrição será | = |x | rá | x, ou —1=-—x. +x, que 
= 6 Ei E 4 2 sa colucou na linha 7 e vem: 
E = 3” X, E x. T X4 OU 
; az ) | Base c E Fa Ea Po Po 
—2=—x—x+x. que se colocou na última do Õ 91 = =] e 0 
quadro anterior. Como a coluna pivô se transforma em 
zeros e as colunas P,, P,, Pg. P', P',, P', não se alte- Ba e! Ez , 18 o 
ram na transformação, podemo-los eliminar dos quadros: 2 P [8] —4 8) mito 1 o) 
3% PB; O 20 3 15 —24 (6) 
i Base c|P, Fi P, F SA 
a e ee 4 PP, 0 3 1 2 == O 
O q = na: A O 
5 Ps, 14 O =| O O O 
RR O| 6 -31 —15 7 O) 
8 PF 4 1 O = 1 O 
=X R O| O O 1 3 O 
7 == O me À 6) 1 
ad o P. 0-4 3] —24 15 Ô 
a equação geradora é a segunda, logo 
4 P,; 410] O 1 =:3 2 8) 
So Pg md 9 + 0 O 0 ii qui — pe) =1 > P, coluna pivô. m,.< m, < 
o Pu 21] & (é) 7 = 6) 
A a a 3 3 
1 1 0 4 0 1 “—s =15A pera = 3. logo a restrição será — 1 
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— 
E 


— X, É Xyy já colocada na linha 7. Transformando vém: 


i Base c Ko P, E, Po 
O 52 =—4 = = 
TP; DB 14 -3 —415 7 
à Pp. 6 2 0 1 —3 
q E 5 3 —24 15 
+ Pt 1 1 = 2 
5 PP ma 0 —1 O) 0 
» PB 2 2 0 1 =1 


+ 


e por não haver qualquer termo independente (da co- 
luna Es ) negativo, chegou-se à solução final inteira 
que é: 

valor óptimo (mínimo) e inteiro 52 


==]; m=x,=0n=x,=2;x=14; 


x, =2;x=5 


6.6 — Programação do Algoritmo de Gomory 


Para facilidade de utilização convém que os dados 
a fornecer ao programa sejam fáceis de obter a partir 
do problema apresentado. Assim, além da programação 
do algoritmo de Gomo:y propriamente dito, desenvol- 
vemos uma função «PRONTA» capaz de submeter a 
mat:iz À das restrições e o vector B dos coeficientes 
da função objectivo ao tratamento necessário para con- 
seguirmos o quadro inicial exigido pelo algoritmo de 
Gomory. Em seguida a função «OPTIMAs, que interpre- 
ta o algo.itmo de Gomory, dá-nos a solução do proble- 
ma ou a certeza da sua irresolubilidade pelo referido 
método, mas por uma questão de exposição e até de 
memorização, compomos estas duas funções numa outra 
de nome «GOMORY» que só por si efectua todos os 
cálculos necessários. 

Se o problema for o exemplo anterior os dados a 
fornecer se.ão apenas 


A—B MIS —1 O O 12 
p) 2 7 O —1 O 14 
9 8 3 O O “1 W 
oO OO 


B<— 10 14 29 


e de forma semelhante se procederia para qualquer pro- 
blema depois de colocado na forma canónica. 

Vejamos então a função «PRONTA» que prepara 
os dados para a função «OPTIMA» e com a explicação 
das variáveis e dos cálculos efectuados em cada linha: 


à A PRONTA B 


Fixação da dimensão da matriz A; multiplicação da 
matriz À por —1 se não tiver elementos independentes 


TÉCNICA 435 


negativos (agora na primeira coluna) e que passa a 
ser AR 


[1] ARSQ(TISI(v/AD1MI DO) XA-IGSA; 
D + (pA) [1]; 6 <— (pA) [2] 


Selecção dos vectores unitários de AR, V, sua elimi- 
nação e definição da nova matriz AC 


[21] AC [((0,9V) < Afe=Vve-t=++41A[;1+ 
+ sG — 1])/AR 
Definição de E, número de colunas 


[Ma Eett-Fe+/V)=1 
B1, coeficientes da função objectivo ordenados por vec- 
tores coluna e com os zeros necessários. 


[ 4] B1< ((—V)/B), (V/B), (D— F)pO 
B3, índices dos vectores da base 


(6) Bae-(E+s Ds É 


Posição dos vecto;es unitários, U; há vectores unitários 
na matriz dada? Salto para a linha [8]. 


[6] =>(1—v/V)/8;U <-(0,V)/s6G 

indices da base para o caso de haver vectores unitários. 

[ 7] B3[spL] <L;B3 [B3s5L] <— B3 [splLe(V/ 
/[s6-1)+.xX — AR [;U]] 


Coeficientes Cc B2, 
da base 


co respondentes aos vectores 


[ 8] B2<-B1[B3] 

Matriz W destinada a ser usada por «OPTIMAs. 

[9] We qWIDp1),EpO)/9 AC) + (ID p O), 
EplN— ADE) =sE 


Matriz auxiliar, argumento de «OPTIMA» composta por: 
coeficientes da função objectivo; coeficientes da base: 
indices dos vecto;es da base. 


[10] Y< (3,0 B1)p B1, B2, B3 
Y 


Vejamos agora a função «OPTIMA» 
7 À OPTIMA B 


Número de vectores da base dual, CD 


[7] CDe-)-k=ipÃ 


Cálculo do vector ZE 6 ZMC 


a 
[2] ZME< (0, — B[llys (p A) [2] — 1]) + 
+ (B[2;]) [sCD]+. X A[:CD;] 
Ampliação da matriz A com a primeira linha z, — c, 
[3] AL] -A[1;] +ZMC;A< (O, (pA)[1] 
p= A 


Há elementos negativos na primeira coluna? Se não sai 
a solução. 


33 


[4 =>(A/0<XEGS-,A[ +slpA) VA PRONTA H 


Ae Latvia D;I]r0)jnÃo I0A De teA)LI Tee inA ILS | 


[13 
LE ACO VIsa/=-Pel=+2]ALLeatoA)LZ)-1])/AR 
= 15: 1]1/14 [3]  Es+(G-Pe+/V)-1 
N F Rar a [5] R3+(E+1D),1E 
Selecção da equação geradora, índice L, da restrição LE] +(L=v/V)/OGUeto,V)/6 
[7] Ralip4J+L;83[831,5])+eRalipie0F/10-1)+,=m=ARC GUI 
de corte, [a]  B2+B1[83] 
pa) WeB(((DEI) EDOJNDAE CL DOO),BoLIA-LO,1E Jo, =1E 
D) Si ; 2,0: 
[6] L&A1+HEs [E;OCATI;]<O; a AR can 
G— (Ne ,E),1 
Ta DPTINA & 
[13 EheL+-/pA 
| PES lacei | [2]  aMge(0,-BLLGA(DAJC2]-1])+(8027])010D]+,xAlaCD:I 
ja elementos negativos na Aqueça seleccionada de E PA nn dd da 
baixo de elementos negativos da linha z; — c, ? Se [4 jo +(a/05E+,ADIeitpAJLIT-191])/14 
Ee - E [5] L+1+EiL/EGQ+ALI;]<0;6+(N+pE),1 
não há não tem solução, [6] +[0=+/Qxe+D, (ALL;])Li+itoM)L2]-1]<0)/17 
[71 PIvOsG,ALL SID MINeL/-FelPeo)/, ALI; 
[8] LOM+*[/U=-(PuqQ)/ALL:])+(L1000+Vi=MIN)-1000 
[6] =>(0=+/0 XP O, (A [L:]) [1+ [9]  LDA+LDA+I=SLDA 
[410] Fi O o io 
ae Lã1il] AstloAIL2JoL)OM 
; : lp A] [2] 1] « 0)/17 [12] ArA+C,ALEIVO) o, “RESCO 
[13] +u,pAc(Dz+dAxA N/A 
| | 3 mivi [44] 'FALOR DPTIMO = !4A[1;1]) 
Determinação da coluna pivô, PIVO. o) Pants = VIRADA ad 
[16] + 
[7 PiIvVO— GALN])s—-TMIN-— /—-Ve- (PX ro FNAO TEN SILUCAO! 
| ci à 
xXx 0)/,A [1:] 
erro do Interpretador na função chão ( ), LDA. Ea) Dad 
Determinação de À com um artifício para contornar um se cs 
[ 8] LDA T/(- (PXO)/A [L;]) + (L 1000 + pesso 
ql. 6] 10 24 24 000 
+ V => — MIN) — 1000 P+3 700 
Plin Tl 6 3920014 
| ; — PLZ 7) 2270101 
[9] LDA < LDA + 1 LDA Plãji 7)+ 963400140 
A JCARD END 
Restrição de corte, RESCO, deduzida de A [L;]. 
[10] RESCO<((L 1000 + A [L;] > LDA) — 1000), 1 
Ampliação da matriz À com uma coluna de zeros. 4: RPE & 
[11] A ((lpA) [2]21),0) A 
Aplicação das fórmulas de eliminação completa, PRONTALI) das 
“12 .B 11 d .g E 
[12] Ae A+ (AL PIVO])o. X RESCO Ss e GQ ray 
Eliminação das colunas de zeros e regresso a [4]. PRONOARES, ss eg 
“1u E “a E. 
Ra A mm 
Saida dos resultados PRONTAL4] 10 14 21 0 0 8 
PRONTALS] 4 5 6 1 2/3 
é mm PRONTALG] 5 6 T 
[14] “VALOR OPTIMO = “;A [151] ONEREA E GM ME 
PRONTA] O O O io th 21 
[15] 'VARIÁVEIS = E [B[3;] ss N] PRENGALAA 4 rg 
Au E. a 7 
[16] +=0 a =] : e 
á D “a - 
A A À ' ú 4 Ú 1 
[17] NÃO TEM SOLUÇÃO sRONTÁLIO) 
io 14 al Q 0 a 
V õ É E 10 Ly 21 
4 5 [7] 1 Ê 4 
E finalmente a função «GOMORY» 
A GOMORY B N 
JE. | “so “pg Cj “a 
[1] A PRONTA B sd (Om MR ug 
“40 9 "o Ts 
[ 2] W OPTIMA Y : : = E 
v | E) Q "1 
8 
Apresentamos em seguida seguida a listagem das 1 q 2 D 5d 
funções pelo computador e, para melhor compreensão ; 
dos algoritmos, apresentamos também a sequência dos 
cálculos efectuados na resolução do problema já resol- 10 14 21 0 Do 
z [i! ] [É 10 lu 1 
vido manualmente. 4 : 6 1 2 3 
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TSOPTINA+118 
| OPTINA F 


OPTINALI) 3 

OPTINA[2] 0 do 

OPTINAL3] O | 10 
0 


DPTINAL4] 
OPTINALS] 30 1 1 15 1 
DPTINAL6]) 
OPTINALT] 2 
DPTINA(8] 3. 
DPTINALI] 9. 
OPTINA[10) 
OPTINAL11] 
a) | 


de om um 


| 
“sá 
14 
10 
[, | Ô Ú 
Í 1 Q 
[| a | 1 
DPFINMA[12] 
Ce] [| 
20 ) 
“6 3 
26 E 
o) u 
[| Ú 
q 0 
OPTINAL13] u 
OPTIMAL4]) 
DETIMALS] 30 1 1 LE L 
DPTIMAL6] 
OPTINALTI) à 
DETINALB] 3 


a ta Em 


1 


oDadDdona 


[O O E 


Oordtoo 


1] Tom 1 1 1 
OPTINA J 
uo ta. 
z0 3 
“6 0 
26 3 1 
! 
1 
õ 
1 
+ 
1 


O RS UM qd dl dé 


ED R-D-D=Denoa, 


DPTIMAL1I3]) dq 
OPTINAL4] 
DPTINALS) 40 


te 
a 
= 
e 
Es 


Capítulo 7 


VARIÁVEIS BOOLEANAS E ALGORITMO «BALAS 


7.1 —— Transformação de Variáveis Inteiras em Variáveis 
Booleanas. 


Há muitos problemas de P.,L.l. cujas variáveis são 
por natureza booleanas mas em muitos outros as va- 
riáveis tomam quaisquer valores inteiros positivos em- 
bora geralmente não sejam muito grandes. Por outro 
lado alguns algorimos trabalham apenas com variáveis 
booleanas e outros embora aceitem qualquer tipo de va- 
riáveis são frequentemente mais eficazes quando são 
booleanas as variáveis envolvidas. 


Ora estes factos fazem com que seja de grande 
utilidade saber transforma; variáveis inteiras quaisquer 
em vailáveis booleanas até po-que com esta transfor- 
mação aumentamos a capacidade do algoritmo que va- 
mos tratar e que passa a poder resolver inclusivamente 
os problemas resolúveis pelo Método de Gomory. O 
único obstáculo que pode su:gir é aparecer com esta 
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DETIMAÇ6E] 

0PTIMALTI] 3 

OPrINAÇO)] 24 

OPTIMAL9] 24 

OPTINAÇ10] "1 O “1 0/4 


DPTINA[11] 
uz Ts “8 | [a] 
6 “a CM 7 e 
o o 1 “3 [E 
“ug 3 24 415 o 
| 1 “3 7 o 
[1 “4 q [o] o 
2 o 1 “4 [o 
OPTIMAL12] ” a 
51 lá [5] 4 E] 
21 “a [ r “45 
= U D CH = 
20 a o] 15 Zu 
3 1 [+] 7 a 
e] E õ 0 Õ 
1 0 0 "4 i 
DPTINÁALIA] 4 7 4 
BPTINA[4] 
DPrINAIS] 30 1 1 16 1 
DPTINAÇE ] 


DPTINAL7] 3 
0PTINALO] 3 
DPTINAL9]) 3 
DPTINA[10] “1 q 1 0 1 
OPTIMAL11] 


51 lá i E] D 
21 3 7 "15 fr 
: ] | 1 E] 
20 | 15 “24 a 
3 af 2 “a O 
o | a] 0 o 
i o “4 1 q 
OPPIMAL12] 
42 “4 0 “8 4 
mm “3 o “415 7 
2 o o 1 3 
5 3 o “zu 15 
1 1 0 “a 2 
0 “4 E] 0 0 
2 õ õ 1 E À 
DPTIMAÇ13] 4 7 4 


DPFINALM) 14 

OPTIMALi4] VALOR OPTIMO = 52 

DPTINAL15] FARIAVEIS = 12 0 27 24 2 5 
OPTINA[16] O 


P GONORF Q 


FALLON DETINO =» 52 
VARIAVEIS = 1 0 2 44 2/5 


trensformação um g'ande número de variáveis boolea- 
nas. Vejamos então como é feita tal transformação: 
Qualquer variável inteira não negativa x se pode 
considerar como pertencente a um certo intervalo 
(O <x<L) onde L é o seu limite superior. E sendo 
p o inteiro positivo tal que 2P”I <L < 2P sabemos que 
qualque: núme.o inteiro se pode escrever sob a forma 
duma soma de potências de dois da seguinte maneira: 


(3.1) 


onde 5, são variáveis booleanas. Temos então a variá- 
vel x substituida por p variáveis booleanas. 

Exemplo: suponhamos que a variável x inteira e 
não negativa toma no máximo o valo; 20, isto é, 


(O < x < 20) com 2º < 20 < 25. Então será: 


H 
x= 4 995- 5; 
1 


+205,=165,+85,1+48,+28,+3, 


=243,+2:8,+ 228,+218,+ 
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soma que dá 0, se 5, =... =5,=0; dá 1.se 9,=1 
e5,=...=65,=0; dá 2 se apenas 5,=1; dá 3 se 
apenas 59,=60,=1 etc. até 31 quando 3,=... = 
=8,=1, sendo 31=2º—1. 


A equação (3.1) define então uma variável que pode 
tomar qualquer valor inteiro entre O e 9p—1, portanto va- 
lores que podem ser até malores que o limite superior 
L de x. Isto faz com que certas combinações de ô, = 1 
se possam por vezes desprezar conforme a posição de 
L relativamente às potências de dois e como se vê nos 
três casos seguintes: 


1.º L=2P —1 (caso de L=31 no exemplo): 
são possíveis todas as combinações de 5 — 1. 


2.º L=2P71 (caso de L= 16): podemos es- 
crever a restrição de implicação: 5, =1 5 6, =0 
para 2<i<p ou, doutraforma: ô,+ó, <7 paa 


2<is<p. 


3.0 2P71 <«L<2P — 1: Começamos por defini: o 
número L' por: 


p= 


Ee=L = :9P 7 t=sO 
donde por serL< 2” vem L' < 2?! Definamos agora 
o inteiro positivo q de modo que 


9p—1 ELE 24 


e portanto q< p— 1 por ser L'< 271, Chamando | 
ao conjunto dos índices i tais que q<p-—iou 
(<isp-— q, podermos escreve; 


L'<2P-1 para qualqueri<p-—q,e portanto 


L<2P-i + 2P-1 para todo i<p-a. 


Então olhando a definição (3.1) 


p—da ; j' ) 
x= 9p-—1 Em ” 9 p-i 2) ap-1 
|=2 i=p—qr 


vê-se que para este limite superior de L, salvo se 
q=p-— 1,se pode escrever a restrição 0, = 1 > 8, = 
=0 ou d,+ô,<1paraZ<si<p-—a. 


E três situações podem oco.rer: 
a) L'=29 — 1, Então a única restrição possível 
á a anterior pois para além disso todas as combinações 


são possíveis, 


b) L'=25 1=2p-fp—9t1), pela definição de 
L*, vem 


L=9P-1 4 9p-p-q+1) 


8 comparando com o segundo caso, vêm as seguintes 
restrições a acrescentar à anterior: 5, = 108 


pad 
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=15> 8,=0 ou 5, +56, + 8 <2 para 


tra | 
p-9+24a sp. 


c) 2P1 <L' <29 1, Define-se o número L"= 
=L'— 297150 e seguem-se os cálculos feitos para 
N' e assim sucessivamente até nos encontrarmos na 
alinea a) ou b) que faz terminar o processo, 

Para ilustrar o que acabamos de expôr, tomemos 
como exemplo o problema pelo algoritmo «Gomo:y» mas 
na seguinte forma: 

Minimize 10x, + 4x, + 21x, sujeito a x < O 
e inteiras, com 

= =17%— 3x E 
- 2% = Emma eta 
-u-— kr 3% &— 10 


Em primeiro lugar temos de encontrar um limite 
supe ior para as variáveis. Embora se veja que esse li- 
mito podia ser 7, tomemos L = 12 para discutimos o 
3.º caso: 


Ogh SIL p=4 x= 
4 

= > 23, =88,+45,+25,+54 
j=1 


Temos o 3.º caso pois 241 <12< 24-71; L'= 
=12-2H1=4, 2HigACZ: q=3, e como 
q=p-— 1 não é possível a primeira restrição do n.º 3. 
Mas L'= 4 = 2º-1, alínea b) e portanto 


d,+t59,+8, <2 para 3<i<4, ou 


à rosto + do, GS 
o que é verdade pois bastam duas variáveis dó =1 


para dar o maior valo: possivel (12) para as variáveis. 
Assim o problema transforma-se em: 


Min 805, + 405, +205,+ 108,4 

+323,+188,+83,+48,+4 

+ 1683, +845,, +428,,+215,, 
— 645,—328,-— 168, -84, - 883, -- 


— 4483,—228,—-118,— 728, — 


Cio E DE E — 12 
STO Toro, < 2 
O Top Tr OS Td; < 2 
do do td Fa s 2 


e .esolvido o problema com as doze variáveis booleanas, 
para obtemos o resultado correcto x =1; x,=0; 
xy = 2 e valo: óptimo 52, teria de ser 3,=8,=1e 
as restantes nulas, ou seja (0, 0,0,1,0,0,0,0,0,0, 
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1, O,) que daria para valor óptimo 10 X 1+42X1= 
= 52, 


7.2. Métodos de «Enumeração Implicitar 


Um problema de P. L. |. de variáveis bcoleanas 
pode ser assim apresentado na sua forma canónica: 


n tm 
min Z = > c; x; com ba a 
j=1 =1 


E | 


41; X ; inteiro na hipótese de c > O, para 


EA Lai 6 E= 12 ces 


e isto porque basta fazer a transformação x = 1-—-x 
para toda a variável x, para transformar a maximização 
em minimização. E se fizermos a mesma transformação 
às variáveis que na função objectivo a minimizar te- 
nham coeficientes negativos, conseguimos que todos os 
seus coeficientes sejam não negativos. 


Chamamos variáveis de folga às variáveis x > 0 
tais que 


n 

W a... x. É =| 
ca É x +x =b 
i= 


e solução ao conjunto (x; 634) de elementos que 
satisfaçam a igualdade anterior, que sejam inteiros e 
tais que pertençam zo intervalo [0; 1]. Existem por- 


tanto 2” soluções. Chama-se solução admissível à 


solução anterior que satisfaça também xt = O e será 
óptima se ainda minimizar a função objectivo. Vejamos 
ago.a em que consiste o princípio de Enumeração Im- 
plicita: 

Caracterizando a solução S, pelo conjunto J, = 
= [1,2,...,n) constituído pelos índices das variáveis 
X; = 1 nessa solução e sendo S, uma otura solução 
tal que J, c J,. chamar-seá a S, solução descen- 
dente de 8, . O princípio dos métodos de enumeração 
implícita é não fazer cálculos senão num sub-conjunto 


das 2" soluções. Assim do conjunto das soluções a 
examinar desprezam-se todas as já examinadas e aque- 
las que embora não examinadas são descendentes dou- 
tras já examinadas, ou por termos a certeza de condu- 
zirem a valo-es da função objectvio superiores ou por 
serem não admissíveis, Do facto de ao ser analizada 
uma solução outras o serem implicitamente, vem o nome 
do Método de «Enumeração Implicita». 


Lego se vê que o melhor método será aquele que 
permita desprezar o maior número de soluções descen- 
dentes, 


Os diversos algoritmos que seguem este princípio 
dife:-em entre si pela maneira de ver se uma determi- 
nada solução já foi examinada, pela escolha da variável 
igual a um para passar a uma solução descendente 
ainda não examinada e pelo processo de voltar atrás 
pa;ra escolher um novo conjunto Jp a examinar. 
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7.3. Algoritmo de Balas 


7.3.1. Generalidades e Notações 


Foi Balas quem apresentou o primeiro algoritmo 
baseado nos métodos de enumeração implícita, Poste- 
riormente apareceram algumas ideias para melhorá-lo 
sob certos aspectos, especialmente por F. Glover com as 
suas restrições especiais ou a chamada «técnica do fil- 
tros e ainda por A. Geoffrion. 

Inicialmente parte-se duma solução S, em que 
o correspondente conjunto dos índices das variáteis 
iguais a um, Jo, » é vazio é vai-se passando a novas 
soluções (ou a novos caminhos possíveis do grafo apre- 
sentado) escolhidas sob certas regras. E assim que se 
atingem certas condições abandona-se o caminho se- 
guido e metemos por outro até à solução óptima no 
caso da sua existência. Antes de entrar propriamente 
no algoritmo convém referir as notações seguintes: 


[Ja b 


«f).3 


+/J]20 e[Jzl + []=2 


onde +/J designa o número de variáveis, ao lado de 
cada ponto, com o valor 1. 


— Às soluções são numeradas por S,, 8,, ..., Sp 
à medida que se vão obtendo e conservam o seu nú- 
me:o ao longo dos cálculos. 


— Designa-se por Z, o valor da função objectivo 
correspondente a S, , com Jp a designar o conjunto 
dos índices das variáveis iguais a 1. 


— Z, ou simplesmente Z será o valor da função 
objectivo correspondente à solução admissível S, . Con- 
venciona-se Z, =+ 0 se não fôr admissível a solução 
encontrada. 


— As variáveis de folga cor-espondentes à solução 


S, serão designadas po; E 


” ou apenas Xp : 


7.3.2. Algoritmo 

1.º Iniciam-se os cálculos tomando para primeira 
solução S, em que nenhuma das variáveis é igual a 
um, ou 8, =10, 0, ...., O). Se após examinada esta 
solução se concluir que todas as outras estão implicita- 
mente examinadas, para-se, se não, passa-se à solução 
S, acrescentando uma variável o RR A paragem é 
decidida pela aplicação da «regra de paragem» do nú- 
mero 3 e a nova variável é escolhida segundo a regra 
da escolha do número 5. 
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2.º Conhecida uma solução terá de aplicar-se a 
«rega de paragem» e para isso vamos definir o conjun- 
to N, dos índices não pertencentes a Jd, (J, é vazio), 
mas que juntos aos deste podem definir «melhores» 
soluções que as anteriores. Para isso começamos por 
considerar todos os índices possíveis (das variáveis 
booleanas) e excluimos os que pertencerem aos se- 
guintes conjuntos: 


— D, constituído por todos os índices | que 
unidos aos de )y dão uma solução cujo correspondente 
valor da função objectivo é maior que Z, (só se con- 
sidera D, se Sp é admissível, para termos ê, finito), 
ou seja: 


D.=Ij:Z 


p ui E; > Z,| 


(atendendo a quo usamos variáveis booleanas). 


— E. , Constituído por todos os indices j que fa- 
cam com que as variáveis de folga negativas em Sp 
continuem ainda negativas mas com maio:es valores 
absolutos nas soluções descendentes de Sp , OU Seja: 


E, = () 18) > O para todo o | dE LAR O; 
pois x! = -a, 
p+1 Lp À.) 
— p' constituído po: todos os indices que unidos 


ao de J, dão origem a uma solução já examinada an- 
tes de 2» - OU seja: 


K p = (| “Ju (j] já foi examinada] 


36 Regra de paragem: Conhecido o conjunto N, 
para-se em S, | isto é, não se examinam as soluções 


descendentes de 9 | se: 


a) Sp é solução admissível (as variáveis de folga 
são não negativas); 


b) N, é vazio, não existindo portanto qualquer 
indice que unido a J » possa dar melhor solução; 


c) No não é vazio mas para uma linha i tal que 
E e - 
Xi p S O tem-se 
” Wo 
K IE Np , 


já que mesmo tornando iguais a zero as variáveis com 
coeficientes a, | > O e iguais a um as variáveis com 
coeficientes à; | < O, a variável de folga continua 
negativa. 

Se não ocorrer nenhum destes três casos passamos 
ao número 5.º para continuar os cálculos com nova so- 
lução descendente de Sp : se algum deles ocorrer pas- 
samos ao número seguinte. 


4.º Regra do Passo Atrás: Consideramos todas 
as soluções do mesmo caminho percoriido de S, as p 
e numeramo-las por ordem inversa: S, o (= Sp). Sa! À 


E a pone |=8, ). Destas tomamos a primeira 


dé ap ) 
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solução Sal cujo correspondente Nas contenha pelo 
menos dois índices (se contiver apenas um não vale 
a pena examinar as descendentes pois foi delas 
que acabámos de sair, conforme se vê no grafo) e de- 


finimos um novo conjunto Na resultante de Na ao 
eliminar-se-lhe: 

a) todos os indices que unidos a dat dão uma 
solução já atingida. 

b) todos os indices que unidos a dat dão uma 


solução cujo valor da função objectivo é superior a Z 
Conhecido este novo conjunio N voltamos a apli- 
car a regra de paragem do número 3.º e, caso não se 
pare, o que obrigaria a recuar ainda mais, passamos ao 
número 5.º para escolher a nova solução descente. 


5. Regra de Escolha: Chegados à solução Sp 
passa-se a uma solução descendente juntando a Jp 
um índice j pertencente ao conjunto N,. Vejamos co- 
mo se escolhe esse índice: 


Para todas as xÉ < O tem-se forçosamente 


bd | 
= min (0;a,,)>0 
JENp é 


a 


Então existindo um EN, tal que 


W | 
xo — = min(0;a ,)-la,,I<Ocoma,, >O 
ip ! s4 il L,] 
JE Np 
só se conseguem x > 0 se x = O e portanto evita- 


-so a escolha do índice |. Ora, designando por r o es- 
colhido, as variáveis de folga tomam os nossos valores 


Assim se considerarmos para cada EN, a quantidade 


mm 


Ra pe E : 
Q = ba) min (xi, — a; :0) 
= 


e se não se anular para nenhum ieN,: a solução 9641 
não pode ser admissível, Escolheremos uma solução o 
mais próxima possível duma solução admissível, isto é, 
em que o valor absoluto das somas das variáveis nega- 
tivas seja o menor possível. Para isso seleccionamos o 
indice r tal que 


mi 

“º E | á 

ai | po min x; p a,| :0) 
ND |. i- 


Esta regra faz com que, se a expressão Q se anular, 
seja escolhido o índice r que torne S +, admissível. Se 
aparecerem vários índices para r será seleccionado o 
que provocar menor valor para a função objectivo. Re- 
sumindo, a regra de escolha consiste em: 


— Se existir solução admissível é escolhido o r 
que dá tal solução, 
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